arXiv:1501.05373v3 [hep-lat] 26 Apr 2015 


FERMILAB-PUB-15-019-T 


B —¥ nlv and B s —>■ Klv form factors and \V u b\ from 2+1-flavor lattice QCD with 
domain-wall light quarks and relativistic heavy quarks 

J. M. Flynn, 1 T. Izubuchi, 2,3 T. Kawanai, 2,3, * C. Lehner, 3 A. Soni, 3 R. S. Van de Water, 4 and O. Witzel 5 ’^ 

(RBC and UKQCD Collaborations) 

1 School of Physics and Astronomy, University of Southampton, Southampton S017 1BJ, UK 
2 RIKEN-BNL Research Center, Brookhaven National Laboratory, Upton, NY 11973, USA 
3 Physics Department, Brookhaven National Laboratory, Upton, NY 11973, USA 
4 Theoretical Physics Department, Fermi National Accelerator Laboratory, Batavia, IL 60510, USA 
5 Center for Computational Science, Boston University, Boston, MA 02215, USA 

(Dated: April 28, 2015) 

We calculate the form factors for B —> nlv and B s —> Klv decay in dynamical lattice Quantum 
Chromodynamics (QCD) using domain-wall light quarks and relativistic b quarks. We use the 
(2+l)-flavor gauge-held ensembles generated by the RBC and UKQCD collaborations with the 
domain-wall fermion action and Iwasaki gauge action. For the b quarks we use the anisotropic 
clover action with a relativistic heavy-quark interpretation. We analyze data at two lattice spacings 
of a m 0.11,0.086 fm with unitary pion masses as light as M n « 290 MeV. We simultaneously 
extrapolate our numerical results to the physical light-quark masses and to the continuum and 
interpolate in the pion/kaon energy using SU(2) “hard-pion” chiral perturbation theory for heavy- 
light meson form factors. We provide complete systematic error budgets for the vector and scalar 
form factors f+(q 2 ) and /o(+) for both B —» nlv and B s —> Klv at three momenta that span the q 2 
range accessible in our numerical simulations. Next we extrapolate these results to q 2 = 0 using a 
model-independent z-parametrization based on analyticity and unitarity. We present our final results 
for f+(q 2 ) and fo(q 2 ) as the coefficients of the series in z and the matrix of correlations between 
them; this provides a parametrization of the form factors valid over the entire allowed kinematic 
range. Our results agree with other three-flavor lattice-QCD determinations using staggered light 
quarks, and have comparable precision, thereby providing important independent cross-checks. Both 
B — > nlv and B 3 —> Klv decays enable determinations of the Cabibbo-Kobayashi-Maskawa matrix 
element \V U + To illustrate this, we perform a combined z-fit of our numerical B —> nlv form-factor 
data with the experimental measurements of the branching fraction from BaBar and Belle leaving 
the relative normalization as a free parameter; we obtain \V u b\ = 3.61(32) x 10 -3 , where the error 
includes statistical and all systematic uncertainties. The same approach can be applied to the decay 
B s —> Klv to provide an alternative determination of \V u b\ once the process has been measured 
experimentally. Finally, in anticipation of future experimental measurements, we make predictions 
for B —» nlv and B a —t Klv differential branching fractions and forward-backward asymmetries in 
the Standard Model. 

PACS numbers: 11.15.Ha 12.38.Gc 13.20.He 14.40.Nd 

I. INTRODUCTION 

Semileptonic S-meson decays play an important role 
in the search for new physics in the quark-flavor sec¬ 
tor. Tree-level decays that occur via charged W-boson 
exchange are used to obtain the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements \V u b\ and \V c b\, while 
flavor-changing neutral-current decays provide sensitive 
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probes for heavy new particles that may enter virtual 
loops. Decays involving r leptons are especially sensitive 
to charged Higgs bosons that arise in many new-physics 
models (see e.g. Ref. [1] and references therein). 

The decays B —> nlv and B s —> Klv probe the quark- 
flavor-changing transition b —► u. In the Standard Model, 
the differential decay rate for these processes in the B( s y 
meson rest frame is given by 
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where P denotes the light pseudoscalar pion or kaon and 
q = (pb —pp) is the momentum transferred to the outgo¬ 
ing charged-lepton-neutrino pair. The vector and scalar 
form factors f + (q 2 ) and fo(q 2 ) parametrize the hadronic 
contributions to the electroweak decay and must be cal¬ 
culated nonperturbatively, such as with lattice QCD. 
Given an experimental measurement of the branching 
fraction and a theoretical calculation of the form fac¬ 
tor (s), these decays enable a determination of the CKM 
matrix element \V ub \. (The contribution from fo^(q 2 ) in 
Eq. (1) can be neglected for light leptons £ = e, p, given 
the current experimental and theoretical precision.) To 
date, both the BaBar and Belle experiments have mea¬ 
sured B(B — > n£v) [2-5], and the experimental uncer¬ 
tainty will continue to improve with the collection of data 
at Belle II. The decay B s —> K£v has not yet been mea¬ 
sured, but we anticipate a result from LHCb in the next 
few years. 

The CKM matrix element V ub places a constraint on 
the apex of the CKM unitarity triangle [6-8]. Its value, 
however, is under scrutiny because of the long-standing 
~ 3 a disagreement between V ub obtained from exclu¬ 
sive B —>• nip decay and \V u b\ obtained from inclusive 
B —» X u tv decays, where X u denotes all charmless final 
states with up quarks [6-11]. The value of V ub can also 
in principle be obtained from leptonic B —* tv decay, 
but the current determination from this process lies in 
between those from exclusive and inclusive semileptonic 
decays, and is not as precise [11]. Further, B —> tv is 
sensitive to charged-Higgs boson exchange, and therefore 
does not provide a clean Standard-Model determination 
of |I4 6 |. Thus the decay B s —► K£v , once measured 
experimentally, will provide an important new determi¬ 
nation of \V ub \. 

In this paper we present a new calculation of the 
semileptonic form factors for B — > ir£v and B s —> K£v 
in (2+l)-flavor lattice QCD. Preliminary results were 
presented in Refs. [12, 13]. This is the second in a se¬ 
ries of B-meson matrix-element calculations that uses 
the same lattice actions and ensembles, and our anal¬ 
ysis follows a similar approach to our earlier work on 
B-meson decay constants [14]. We use the gauge-field 
ensembles generated by the RBC and UKQCD collabo¬ 
rations with the domain-wall fermion action and Iwasaki 
gluon action which include the effects of dynamical u, d, 
and s quarks [15, 16]. For the bottom quarks, we use the 
Columbia version of the relativistic heavy-quark (RHQ) 
action introduced by Christ, Li, and Lin in Ref. [17], 
with the parameters of the action that were obtained 


nonperturbatively in Ref. [18]. We renormalize the lat¬ 
tice heavy-light vector current using the mostly nonper- 
turbative method introduced in Ref. [19], in which we 
compute the bulk of the matching factor nonperturba¬ 
tively [14, 16], with a small correction, that is close to 
unity, evaluated in lattice perturbation theory [20, 21]. 
We also improve the lattice heavy-light current through 
0(a s a). 

We analyze data on five sea-quark ensembles with uni¬ 
tary pions as light as ~ 290 MeV and two lattice spac- 
ings of a ~ 0.11 and 0.086 fm. We simultaneously 
extrapolate our numerical results to the physical light- 
quark masses and to the continuum and interpolate in 
the pion/kaon energy using SU(2) “hard-pion” chiral 
perturbation theory (%PT) for heavy-light meson form 
factors [22, 23], which applies when the pion/kaon en¬ 
ergy is large compared to its rest mass. For B —> n£v 
( B s —» K£v), we directly simulate in the momentum re¬ 
gion gLx > <7 2 ^ 19.0 GeV 2 (g^ ax > q 2 > 17.6 GeV 2 ). 
Both statistical errors and discretization errors increase 
at lower q 2 , which corresponds to larger pion/kaon en¬ 
ergies. To extend our results beyond the momenta ac¬ 
cessible in our simulations, we extrapolate our results 
to q 2 = 0 using a model-independent z-parametrization 
based on analyticity and unitarity [24, 25]. Our results 
can be combined with current and future experimental 
measurements of the experimentally measured B —> n£v 
and B s —> K£v branching fractions to obtain the CKM 
matrix element \V ub \. 

There are two earlier published (2+l)-flavor calcu¬ 
lations of the B —> iriv semileptonic form factor in 
the literature by the HPQCD and Fermilab/MILC col¬ 
laborations [26, 27]; updates of these works are in 
progress [28, 29]. In addition, HPQCD recently ob¬ 
tained the first results for the B s —> Ktv form factor 
in Ref. [30]. Both groups use the MILC collaboration’s 
asqtad-improved staggered gauge-field ensembles [31, 32], 
so their results are somewhat correlated. The differ¬ 
ences between the two sets of calculations lie in the 
choices of light valence- and &-quark actions. For the 
b quarks, HPQCD uses the NRQCD action [33] while 
Fermilab/MILC uses a relativistic formulation similar to 
ours. Specifically, they use the Fermilab interpretation of 
the isotropic clover action [34] with the tadpole-improved 
tree-level value of the clover coefficient csw- The more 
recent HPQCD calculation uses the HISQ action for the 
light valence quarks to reduce taste-breaking discretiza¬ 
tion effects, while in the other work asqtad valence quarks 
are used. 
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Our form-factor calculation with domain-wall light 
quarks and RHQ b quarks has the advantage that dis¬ 
cretization errors from the light quarks and gluons are 
simpler, such that the SU(2) heavy-light meson xPT 
expressions are continuum-like. Further, as compared 
to the Fermilab/MILC calculation, we tune the coeffi¬ 
cient of the clover term in the 6-quark action nonper- 
turbatively and improve the heavy-light vector current 
through 0(a s a ), whereas Fermilab/MILC only improve 
it through 0(a). Thus, for similar values of the lattice 
spacing, discretization errors from the heavy-quark ac¬ 
tion and current are smaller in our calculation. Our new 
results for the B —> nils and B s —> Klv form factors 
therefore enable important independent determinations 
of the CKM matrix element \V u b\. 

This paper is organized as follows. Section II provides 
an overview of the lattice calculation. First we define the 
needed matrix elements and form factors in Sec. II A. 
Next we present the lattice actions and parameters in 
Sec. IIB. Then, in Sec. IIC we describe the renormaliza¬ 
tion and improvement of the heavy-light vector current 
operator. Section III presents the numerical analysis. 
First, in Secs. Ill A and IIIB we fit lattice two-point 
and three-point correlators to extract the needed me¬ 
son masses and matrix elements, respectively. Then, in 
Sec. Ill C we extrapolate our numerical data to the physi¬ 
cal light-quark masses and continuum, and interpolate in 
the pion/kaon energy, using SU(2) hard-pion xPT. Sec¬ 
tion IV provides complete error budgets for f+{q 2 ) and 
f 0 (q 2 ) at three momentum values that span the range ac¬ 
cessible in our numerical simulations; for clarity, we dis¬ 
cuss each source of systematic uncertainty in a separate 
subsection. In Section V we extrapolate our form-factor 
data to q 2 = 0 using a model-independent z parametriza- 
tion. We present our results for f+(q 2 ) and /o(? 2 ) as the 
coefficients of the series in z and the matrix of correla¬ 
tions between them; this provides a model-independent 
parametrization of the form factors valid over the entire 
allowed kinematic range. We illustrate the phenomeno¬ 
logical utility of our form-factor results in Sec. VI. First, 
in Sec. VIA, we perform a combined 2 -fit of our numer¬ 
ical B — ► 7 xiv form-factor data with the experimental 
measurements of the branching fraction from BaBar and 
Belle to determine |I4&|- Next, in Sec. VIB, we make 
predictions for Standard-Model observables for the de¬ 
cay processes B —► -k£v and B s — ► I\ iv with £ = /j,t in 
anticipation of future experimental measurements. Sec¬ 
tion VII concludes with a comparison of our results with 
other lattice determinations, and with an outlook for the 
future. 


II. LATTICE CALCULATION 


Here we present the setup of our numerical lattice cal¬ 
culation. 


A. Form factors 


The B —► i t£v and B s —>• Kiv semileptonic form factors 
parametrize the hadronic matrix element of the b —> u 
vector current = u^b: 


(P|V A ‘|I? (s) ) = f+(q 2 ) ( Pb { s) + Pp ~ 


M 2 - Ml 




fo(q 2 )- 


M 


B, 


(*) 


— Ml 


( 2 ) 


where f+{q 2 ) and fo{q 2 ) are the vector and scalar form 
factors, respectively. It is convenient in lattice simula¬ 
tions to instead calculate the form factors f\\(Ep) and 
f±(Ep), which are defined by 


<P|V'*|i? w > = ^2 M B(3) [v»f\\{E P ) +p^f ± (E P )\ , (3) 


where E P is the outgoing light pseudoscalar meson en¬ 
ergy, v M = p^^JMb { s) is the P( s )-meson velocity, and 
= Pp — ( p P ■ In the R^-meson rest frame, 

which we will use for our simulations, /|| and f± are pro¬ 
portional to the hadronic matrix elements of the tempo¬ 
ral and spatial vector currents: 


f\\(Ep) = 


(P|V°|B( B) ) 

V 2 Mb (s) 


f±(E P ) 


(P\V l \B {s) ) 1 
^J2Mb m p p 


(4) 

(5) 


The vector and scalar form factors can be easily obtained 
from /|| and f± via 


U(q 2 ) = 


^2M l 


'(«) 


f\\{Ep) + (Mb (<) - Ep)fj_(E P ) 


fo(q 2 ) = 


2 \ _ \/ 2 Mb( s) 


M 2 B(s) -M 2 p I 


(M Bw -E P )f n (E P ) 


+ ( E P — Mp ) f±(E P ) 


( 6 ) 


(7) 


B. Actions and parameters 

We use the (2 + l)-flavor domain-wall fermion and 
Iwasaki gauge-field ensembles generated by the RBC and 
UKQCD collaborations [15, 16]. We perform measure¬ 
ments at five different light sea-quark masses mi and at 
two lattice spacings of a « 0.11 fm (a -1 « 1.729 GeV) 
and a ~ 0.086 fm (a -1 « 2.281 GeV). The light sea-quark 
masses nil correspond to pion masses of 289 MeV < 
Mir ^ 422 MeV. The up and down sea-quark masses are 
degenerate and the strange sea-quark mass is tuned 
within 10% of its physical value. The spatial volumes 
are approximately (2.6 fm) 3 , such that > 4. We 

summarize the simulation parameters in Table I. 
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TABLE I. Lattice simulation parameters [15, 16]. The columns list the lattice volume, approximate lattice spacing, light (mi) 
and strange (m/,) sea-quark masses, residual chiral symmetry breaking parameter m Tes , physical u/d- and s-quark mass, unitary 
pion mass, number of configurations analyzed and number of sources. The tildes over afh u / d and afh s denote that these values 
include the residual quark mass. 


(f) 3 x(?) 

s» a(fm) 

a -1 [GeV] 

ami 

am h 

(27?7-res 

am u / d arris 

M,r [Me V] 

# configs. # time sources 

24 3 x 64 

0.11 

1.729(25) 

0.005 

0.040 

0.003152 

0.00136(4) 0.0379(11) 

329 

1636 1 

24 3 x 64 

0.11 

1.729(25) 

0.010 

0.040 

0.003152 

0.00136(4) 0.0379(11) 

422 

1419 1 

32 3 x 64 

0.086 

2.281(28) 

0.004 

0.030 

0.0006664 

0.00102(5) 0.0280(7) 

289 

628 2 

32 s x 64 

0.086 

2.281(28) 

0.006 

0.030 

0.0006664 

0.00102(5) 0.0280(7) 

345 

889 2 

32 s x 64 

0.086 

2.281(28) 

0.008 

0.030 

0.0006664 

0.00102(5) 0.0280(7) 

394 

544 2 


In the valence sector we use for the light quarks the 
domain-wall action [35, 36] and generate propagators 
with periodic boundary conditions in space and time and 
with the same domain-wall height (M 5 = 1 . 8 ) and extent 
of the fifth dimension ( L s = 16) as in the sea sector. 
We generate both unitary light valence-quark propaga¬ 
tors with the same mass as the light sea quarks and prop¬ 
agators with a mass close to the physical strange quark. 
On the coarser ensembles we choose am s = 0.0343 and 
on the finer ensembles am s = 0.0272. 

For the bottom quarks, we use the Columbia version 
of the relativistic heavy quark (RHQ) action [17] to con¬ 
trol heavy-quark discretization errors introduced by the 
large lattice 6 -quark mass. We use the anisotropic O(a) 
improved Wilson-clover action with the following three 
parameters: the bare-quark mass mod, clover coefficient 
cp, and anisotropy parameter £. In this work we use the 
RHQ parameters tuned nonperturbatively in Ref. [18] to 
reproduce the experimentally measured R s -meson mass 
and hyperfine splitting; we list their values in Table II. 

We reduce autocorrelations between our lattices by 
shifting the gauge fields by a random 4-vector before cre¬ 
ating the sources for the valence-quark propagators used 
in the 2-point and 3-point correlation functions. This 
random 4-vector shift is equivalent to placing the sources 
at random positions in spacetime but simplifies the sub¬ 
sequent analysis. On the finer ensembles, we double the 
statistics by using two sources per configuration sepa¬ 
rated by half the lattice temporal extent. 


C. Operator renormalization and improvement 

To match the lattice amplitudes to the continuum ma¬ 
trix elements, we multiply by the heavy-light renormal¬ 
ization factor Z'S : 

<P|V <1 |B w ) = Z“<P|y M |B ( . ) ), (8) 

where and V), are the continuum and lattice current 
operators, respectively. Following Ref. [19] we calculate 
the renormalization factor Zy using a mostly nonpertur- 
bative method in which we express Zy as the following 


product: 

Zv, = PvJ Z v Zl v- (9) 

Most of the heavy-light current renormalization comes 
from the flavor-conserving factors Zy and Zy. The re¬ 
maining factor py is expected to be close to unity because 
most of the radiative corrections, including contributions 
from tadpole graphs, cancel [37]. 

Both flavor-conserving renormalization factors Zy and 
Zy were computed nonperturbatively in previous works. 
We computed Zy for our earlier calculation of B-meson 
decay constants from the matrix element of the 6 —► 6 
vector current between two B s mesons [14]. We can also 
take advantage of the fact that for domain-wall fermions 
Zy = Z l \ up to corrections of 0(am res ) and use the de¬ 
termination of Z l \ from Ref. [16]. The flavor off-diagonal 
renormalization factor py is calculated at 0(a s ) in mean- 
field improved lattice perturbation theory [38] and eval¬ 
uated at the MS coupling a^ s (/r = a -1 ). Our perturba¬ 
tive computation extends the work of Ref. [39] to bilin- 
ears with one relativistic heavy quark in the Columbia 
formulation and one domain-wall light quark. For a^ s , 
we use Eq. (167) of Ref. [39], which does not take into ac¬ 
count sea-quark effects. Because sea-quark effects enter 
at two loops, however, and the rest of the computation 
is performed at one loop, the error introduced by setting 
Nf = 0 is of the same size as other remaining truncation 
errors. Further details of the perturbative calculation will 
be provided in a forthcoming publication [21]. Table III 
shows the renormalization factors used in this work. 

We reduce discretization errors in the heavy-light vec¬ 
tor current by improving it through 0(a s a). We compute 
the matrix element of the tree-level heavy-light vector 
current 

V°(x) =q{x)^Q{x), (10) 

plus matrix elements of these additional single-derivative 
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TABLE II. Tuned RHQ parameters on the 24 3 and 32 s ensembles [18]. The errors listed for moa, cp, and ( are from left 
to right: statistics, heavy-quark discretization errors, the lattice scale uncertainty, and the uncertainty in the experimental 
measurement of the B s -meson hyperfine splitting, respectively. 


moa cp ( 

a « 0.11 fm 8.45(6)(13)(50)(7) 5.8(1)(4)(4)(2) 3.10(7)(11)(9)(0) 

a « 0.086 fm 3.99(3)(6)(18)(3) 3.57(7)(22)(19)(14) 1.93(4)(7)(3)(0) 


TABLE III. Operator renormalization factors and improvement coefficients. The flavor-conserving Z factors were obtained 
nonperturbatively [14, 16]. The p factors and improvement coefficients c" were computed at one loop in mean-field improved 
lattice perturbation theory and are evaluated at a^ IS (a 1 ) [21], 



ryll 

Zj V 

rybb 

Zj V 

_ MS / —1\ 
Q! s \CL j 

Pv 0 

PVi 

Ct 

Ct 

cl 

c 2 

c 3 

ci 

a « 0.11 fm 
a « 0.086 fm 

0.71689(51) 

0.74469(13) 

10.039(25) 

5.256(8) 

0.23 

0.22 

1.02658 

1.01661 

0.99723 

0.99398 

0.0558 

0.0547 

-0.0099 

-0.0095 

-0.00079 

-0.0012 

0.0018 

0.00047 

0.0485 

0.0480 

-0.0033 

-0.0020 


operators 


v p( x ) = q( x )^pQ( x ): 

(11) 

v 2 (x) = q(x) 2< DpQ( x ), 

(12) 

v pi x ) = q(x)2j IJ ,-/ i 'd i Q(x), 

(13) 

v pi x ) = q{x)2^^*BiQ{x), 

(14) 


where the covariant derivatives are defined by 
^pQ(x) = ^(U l i{x)Q{x + fl) 


- Ufa - fi)Q(x - £)), 

(15) 

\{q{x + (i)Ul(x) 


-q{x - (i)UX - /})). 

(16) 


The temporal and spatial (H(a)-improved vector-current 
operators are given by the following sums: 

K mp (z) = Vq 0 (x) + c 3 X(x) + cX(x), (17) 

Vr P ( x ) = v °( x ) + c l v H x ) + c 2 V 2 ( x ) 

+c 3 X(x) + cX(x). (18) 

We calculate the coefficients c" and c™ at one loop 
in mean-field improved lattice perturbation theory [21]; 
the values of the coefficients evaluated at a^ IS (a _1 ) are 
shown in Table III. 


III. ANALYSIS 

Here we present our determinations of the form factors 
f+{q 2 ) and fo{q 2 ) for B -A n£u ( B s -A K£u) at large 
values of q 2 > 19.0 GeV 2 ( q 2 > 17.6 GeV 2 ) accessible in 
our numerical simulations. 

Our analysis proceeds in three steps: First, in 
Sec. Ill A, we fit the pion, kaon, and -B(,q-meson 2-point 
correlation functions to obtain the ground-state meson 
masses. The results for these meson masses then enter 
our 3-point correlator fits in Sec. IIIB to obtain the lat¬ 
tice form factors f\\(Ep) and f±(Ep ) at fixed values of 
the pion/kaon energy Ep. In Sec. IIIC, we interpolate 
the renormalized values for fu(Ep) and f±(Ep) in en¬ 
ergy, and extrapolate to the physical light-quark masses 
and the continuum limit, using SU(2) hard-pion \PT for¬ 
mulated for heavy-light mesons. To avoid possible biases 
due to analysis choices, we use the same fit functions in 
the correlator and chiral fits for both processes B -a tt£v 
and B s -a K£v, and fitting ranges that are as close as 
possible. 

We propagate statistical errors throughout the analysis 
via a single-elimination jackknife procedure. We avoid a 
direct dependence on the lattice scale by carrying out our 
analysis in units of the H s -meson mass. The B s - meson 
mass plays a special role because we tuned parameters 
of the 6-quark action to match the experimental value. 
Thus we can obtain the form factors in physical units 
after the chiral-continuum extrapolation by multiplying 
by M^ xp ' to the appropriate power. With this approach, 
the uncertainty on the lattice scale enters only indirectly 
via the values of the RHQ parameters. 

We use the Chroma software library for lattice QCD to 
compute our numerical data for the lattice 2-point and 
3-point correlation functions [40]. 
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TABLE IV. Time ranges used in two-point and three-point fits to determine the lattice meson masses and form factors. For 
the three-point fits, we use the same range for all operators and momenta. 


min, Imax 



2-point fits 

Mk Mb M Bs 

3-point fits 

pBtt rBir rB s K nB s K 

J || j- L J\\ J L 

a at 0.11 fm 
a « 0.086 fm 

[12,23] [12,23] [7,30] [10,30] 
[16,30] [16,30] [9,30] [13,30] 

[6,10] [6,10] [6,10] [6,10] 
[8,13] [8,13] [8,13] [8,13] 


A. Two-point correlator fits 


To obtain the lattice By) —► P amplitude, we first 
calculate the following two-point correlation functions: 

Cp(t,pp) = Y / ^ P ' S (°p&t)Op( 0,0)), (19) 

X 

= 'E(O^ 3) {S,t)d B{B) (jS,0)), ( 20 ) 

X 

Cs w {t) = E<^B w ( f ’*)^ w ( 0,0)), (21) 


where Op = q^^q and Ob { b) = Qj 5 q are interpolating 
operators for the light pseudoscalar and .B( s \-meson, re¬ 
spectively. Both pions and kaons are simulated with a 
point source and point sink, whereas 6-quark propagators 
are generated with a gauge-invariant Gaussian smeared 
source [41, 42] to reduce excited state contamination. 
We employ the same smearing parameters optimized in 
Ref. [18] and denote a smeared source in Eqs. (19)-(21) 
with a tilde above the operator. 

We obtain the pion or kaon energy and B( s ymeson 
mass from the exponential decay of the correlators in 
Eqs. (19) and (20). The correlator in Eq. (21) is used to 
normalize the By) —► P three-point function. We work in 
the Byy meson rest frame such that only pions or kaons 
carry nonzero momentum. In our analysis we use data 
with discrete lattice pion momenta through 27 t(1, 1,1 )/L 
and kaon momenta through 27r(2,0,0)/L. We average 
the results for all equivalent momenta, i. e. with different 
spatial directions but the same total \pp\- We effectively 
double our statistics by folding the two-point correlators 
at the temporal midpoint of the lattice, thereby averaging 
forward- and backward-propagating states. 

At sufficiently large lattice times, the ground-state 
masses and energies can be determined from simple two- 
point correlator ratios. We define the light pseudoscalar- 
meson effective energy and Byymeson effective mass as 


Ep(t,pp) = cosh 1 
Mb (t) = cosh” 1 


Cp{t,pp) + Cp(t + 2 ,p P ) 
C P (t + 1 ,pp) 

Gb( 5 ) ( t ) + (t + 2) 


C B {s) {t + 1 ) 


( 22 ) 

(23) 


We perform correlated, constant-in-time, fits to these 
expressions, choosing fit ranges without visible excited- 
state contamination that lead to acceptable p values. Fig¬ 
ure 1 shows example meson-mass determinations on our 


fine ensemble with ami = 0.004. To minimize bias, we 
use the same fit range for all ensembles with the same 
lattice spacing (although different for light and heavy- 
light mesons); these fit ranges are given in Table IV. The 
resulting pion/kaon and Byymeson masses on all ensem¬ 
bles are given in Tables XV and XVI, respectively. 

In the continuum limit, the pion and kaon energies 
should satisfy the dispersion relation Ep = Mj> + pp 
and the amplitudes of the two-point functions Zp = 
|(O|0p|7r)| should be independent of the momentum pp. 
We obtain the amplitudes from correlated plateau fits to 


Z P {t) 


2 EpCp(t,p P ) 

P-Ept 


(24) 


using the same fit ranges as for the masses. Figure 2 com¬ 
pares the measured pion and kaon energies and ampli¬ 
tudes with continuum expectations on the a ~ 0.086 fm, 
ami = 0.004 ensemble. The measured kaon energies 
and amplitudes agree remarkably well with the predic¬ 
tions from the continuum dispersion relation, to within 
5% even at the largest momentum apx = 27r(2, 0,0 )/L. 
Although the pion data is not precise enough to draw 
strong quantitative conclusions, the measured energies 
and amplitudes still agree with continuum expectations 
within the large statistical uncertainties for all momenta. 
Dispersion-relation plots for the other ensembles show 
similar behavior. 

The kaon data, for which both the energies and am¬ 
plitudes are statistically well resolved, provides an ac¬ 
curate measure of momentum-dependent discretization 
effects, while the pion data provides only a rough cross¬ 
check. On all ensembles, the measured pion and kaon 
energies both agree within statistical errors with the pre¬ 
dictions from the continuum dispersion relation, and the 
measured pion and kaon amplitudes agree with the zero- 
momentum result. Thus, in our determinations of the 
lattice form factors /y and /_j_ in the next section, we use 
pion and kaon energies calculated from the continuum 
dispersion relation (rather than the measured values) to 
reduce the statistical uncertainties. Although we do not 
use the amplitudes obtained from Eq. (24) in our subse¬ 
quent form-factor determinations, the observed momen¬ 
tum independence of Zp provides further support for this 
strategy. 
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FIG. 1. Effective masses of the pion (upper left), kaon (bottom left), B meson (upper right) and B a meson (bottom right) on 
the a ta 0.086 fm ensemble with ami = 0.004. Shaded bands show the correlated fit results with jackknife statistical errors over 
the fit ranges used. All results are shown in lattice units. 




{aKf 


FIG. 2. Comparison of pion (top) and kaon (bottom) energies (left) and amplitudes (right) with continuum-limit expectations 
on the a ~ 0.086 fm ensemble with ami = 0.004. The dashed lines show a power-counting estimate of the leading 0((ap) 2 ) 
momentum-dependent discretization errors. 










































FIG. 3. Three-point correlation function used to obtain the 
B —> P form factors. The single and double lines correspond 
to light- and b-quark propagators, respectively. For B —> nhs, 
the spectator-quark mass (mr) is the same as the light sea- 
quark mass, while for B s —> Khs, the spectator-quark mass is 
close to the physical m s . The light daughter-quark mass ( mi ) 
is always equal to to the light sea-quark mass. Black and grey 
circles denote local and smeared operators, respectively. 

B. Three-point correlator fits 

To extract the desired B ( s ) —>• P hadronic amplitudes, 
we calculate the following three point correlation func¬ 
tions: 

C^{t, Unk,P p) 

= E e *' J <^ w ( 2 > t » k) rp(»^) o p( 5 >0))- < 25 ) 

£,y 

where the improved lattice temporal and spatial lattice 
vector currents V)) mp are defined in Eqs. (17) and (18). 
As shown in Fig. 3, we fix the location of the pion or 
kaon at the temporal origin and the location of the B^ 
meson at time f sn k, and vary the location of the current 
operator over all time slices in between. In our calcula¬ 
tions, the mass of the light daughter quark ( l) is always 
equal to the light sea-quark mass. For B —>■ n decay, 
the spectator-quark mass (l') also equals the light sea- 
quark mass. For B s —> I\ decay, the spectator-quark 
mass is close to that of the physical strange quark. We 
use a Gaussian-smeared sequential source for the b quark 
in the !3( s ) meson to reduce excited-state contamination. 
We insert discrete nonzero momentum at the local cur¬ 
rent operator through p v = 27 t( 1 , 1,1 )/L for B —> n and 
Pk = 27 t( 2 , 0 , 0 )/L for B s — ► K (recall that the 
meson is at rest). To improve statistics, we compute 
the three-point correlators with both positive and nega¬ 
tive source-sink separations (±f sn k); we also average over 
equivalent spatial momenta. 

The lattice form factors /y at and are obtained from 
the following ratios of correlation functions far away from 
both the pion/kaon source and the B( s )-meson sink: 

/|j at (£p) = lim R3,o{t,t snk ,p P ), (26) 

f±\Pp) = 0<< 1™ k i-R3,*(Msnk,£p), (27) 



FIG. 4. Unimproved three-point ratio R^o for p n = 0 with 
several source-sink separations t an k on the a ss 0.086 fm en¬ 
semble with ami = 0.004. 

with 

R 3 ,/i(^5 Gnk> Pp) 

_ C™V(t,t snk ,p P ) J 2 Ep 

IrPn - 

yC£(t,p P )C 2 ; (f sn k-t) 

(28) 

where we use the continuum dispersion relation and mea¬ 
sured light pseudoscalar-meson mass Mp to construct the 
energy Ep. To determine the optimal source-sink separa¬ 
tion for tgnk) ftp) j we carried out a dedicated study. 

We computed the unimproved ratio R® „ for several val¬ 
ues of the source-sink separation on one a « 0.086 fm 
and one a ss 0.11 fm ensemble, choosing those with the 
lightest sea-quark mass because they are most sensitive 
to excited-state contamination. Figure 4 shows the ra¬ 
tio i?® 0 for B st 7r with p n = 0 for several source-sink 
separations on the a ss 0.086 fm ensemble. All plateaus 
overlap within statistical uncertainties in the region far 
from both the source and the sink. The results for the 
ratios R 3 0 and i?® i at nonzero momenta and on the 
a « 0.11 fm ensemble look similar. Because the statisti¬ 
cal errors increase with larger source-sink separation, we 
chose t sn k = 26 ( 20 ) on the a ss 0.086 fm (a ~ 0.11 fm) 
ensembles. This corresponds to approximately the same 
physical distance for the two lattice spacings. 

Figure 5 shows the O(a s a)-improved ratios i? 3 : o and 
Rs,i/p l p for different momenta on the a ss 0.086 fm en¬ 
semble with ami = 0.004. Results for other ensem¬ 
bles look similar. We perform correlated, constant-in¬ 
time, fits to these ratios using fit ranges without visible 
excited-state contamination that lead to acceptable p val¬ 
ues. To minimize bias, we use the same fit range for all 
momenta and ensembles with the same lattice spacing; 
these fit ranges are given in Table IV. The complete fit 
results for the three-point ratios are given in Tables XVII 
and XVIII. 

Finally, we obtain the renormalized B( s \ —> Phs form 
factors / 1 | and f± in the continuum after multiplying by 
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FIG. 5. 0(a s a)-improved ratios db.i/pp (left) and Rzp (right) with t sn k = 26 on the a ss 0.086 fm ensemble with ami = 0.004. 
Plots for B —> nil/ are on the top and B a —> Klv are on the bottom. Fit ranges and fit results with jackknife statistical errors 
are shown as horizontal bands. 


the heavy-light renormalization factors Zy given in Ta¬ 
ble III: 

/||(pp) = 4o/|| at (Pp)> ( 29 ) 

f±(p P ) = Zy/f(p P ). (30) 

C. Chiral-continuum extrapolation 

We extrapolate the renormalized lattice form factors 
to the physical light-quark mass, and interpolate in the 
pion or kaon energy using next-to-leading order (NLO) 
SU(2) chiral perturbation theory for heavy-light mesons 
(HMyPT) in the “hard-pion” limit. In the SU(2) the¬ 
ory, the strange-quark mass is integrated out, and only 
the light-quarks’ degrees-of-freedom are included. There¬ 
fore the chiral logarithms for B —> nlv ( B s —> Klv) 
depend on the pion mass and the pion (kaon) energy. 
The SU(2) low-energy constants depend upon the value 


of the strange-quark mass, as well as on the value of 
the 6-quark mass for R-meson form factors. “Hard-pion” 
yPT, which was introduced by Flynn and Sachrajda for 
the light-pseudoscalar-meson decay K —>• nlv in Ref. [43] 
and later extended to heavy-light-meson decays by Bi- 
jnens and Jemos in Ref. [23], applies in the kinematic 
regime where the pion or kaon energy is large compared 
to its rest mass. Almost all of our lattice simulation data 
is in this hard-pion (or kaon) regime. We can obtain the 
expressions for the B —> nlv and B s —> Klv form factors 
in hard-pion/kaon yPT by taking the limit of the con¬ 
tinuum expressions from Ref. [22] as M^/Ep —>• 0, where 
P = 7T, K denotes the final-state pseudoscalar meson. 

The NLO SU(2) yPT full-QCD expressions for the 
B —> nlv and B s -> Klv form factors in the hard- 
pion/kaon limit are functions of the pion mass Af, r , pion 
or kaon energy E P , and lattice spacing a. They have two 
general forms: 























































10 


/, 


00 J 


no pole 


(M n ,E P ,a 2 )=^ 


n P 


1 + 


/pole’ (M„,E P ,a') = 


Ep 


_ C W 

A p 


^ll 

(4tt/) 2 
1 + 


+ c 


( 2 ) 


' n P \ ' T 


np A 2 

^ , c (2) 
(47t/) 2 p 


Ml 

7r 

A 2 


+ c 


.(3) 


n P A 2 
E 


+ c 


(5)__ 

np A 2 C 


'32 


+ c (4) ^P_ 
A p A 2 


+ c; 


(5). 


A 2 < 


'32 


(31) 

(32) 


one with a pole at Ep = —A = Mb* — Mb ( s) and 
one without. Here the B* resonance corresponds to 
a state with flavor bu and quantum numbers J p = 
0+ for / 1 | and 1“ for f±. The experimentally mea¬ 
sured vector-meson mass is Mp* = 5.3252(4) GeV [10]. 
The scalar B* meson has not been observed experimen¬ 
tally, but its value has been estimated theoretically us¬ 
ing heavy-quark and chiral-symmetry arguments to be 
Mp*( 0 + ) = 5.63(4) GeV [44], while the 0 + -0 _ split¬ 
ting has been estimated in (2+l)-flavor lattice QCD to 
be Mb*(0 + ) — Mp ~ 400 MeV [45]. In our chiral- 
continuum extrapolations we include the effects of res¬ 
onances below the Bn and B S K production thresholds, 
i.e. q 2 < {Mp( s) + Mp) 2 . For B —* niv, the B* meson 
lies below the Bn production threshold, so we include a 
pole in the fit for ff* taking = 45.78 MeV from 
experiment [10]. The predicted value of Mp*( 0 + ) is well 
above Mb + M n , however, so we do not include a pole in 
the fit of /® 7r . For B s —> K , both Mp- and Mp-( 0 + ) are 
below Mp s +Mk , so we include a pole in the fits for both 
/® sK and ff aK , taking A^ sK = —41.6 MeV from exper¬ 
iment [10] and taking Aj?”^ = 263 MeV from the model 

estimate in Ref. [44]. The precise value of Mp*( 0 + ) has 
little impact on the fit because the pole location is so far 
outside the semileptonic region, but we vary its value by 
a generous amount when estimating the chiral-continuum 
extrapolation error in Sec. IV A. 

The one-loop chiral logarithms are the same for / 1 | and 
/_l, but differ for B — niv and B s —> Kiv: 

Sf B * = ~\ (3 g\ + 1) Ml log (jp) (33) 

5f B ‘ K =-\Ml\o g(^) , (34) 

where gb is the B*Bn coupling constant. At tree level, 
the mass of a pion composed of two domain-wall quarks 
is given in terms of the light-quark mass by 

Ml = 2 fl{mi + TOres) , (35) 

where /r is a leading-order low-energy constant. 

We include a term proportional to a 2 in the chiral fit 
functions Eqs. (31) and (32) to account for the domi¬ 
nant lattice-spacing dependence. To make the a 2 ana¬ 
lytic term dimensionless with an expected coefficient of 
0(1) in \PT, we normalize it using the lattice spacing 
on the finer 32 3 ensembles 032 - Discretization errors from 
the domain-wall and Iwasaki actions are of O (oAqcd) 2 ; 


TABLE V. Constants used in the chiral and continuum ex¬ 
trapolations of the B —> niv and B s —> Kiv form fac¬ 
tors [10, 16, 51]. 


a : 

w 0.11 fm ss 0.086 fm 

a -1 1.729 GeV 2.281 GeV 

a(i 

2.348 1.826 

u 

130.4 MeV 

9b 

0.57 

A x 

1 GeV 


using Aqcd = 500 MeV, 1 we estimate these to be about 
5% on the 32 3 ensembles. The remaining discretization 
errors - light-quark and gluon discretization errors in the 
heavy-light current, and heavy-quark discretization er¬ 
rors from both the action and current - are expected 
from power counting to be much smaller. In Secs. IV E 
and IVF, we estimate their sizes to be below 2%. We 
therefore expect light-quark and gluon discretization er¬ 
rors from the action to dominate the scaling behavior of 
the form factors, such that including an a 2 term in the fit 
will largely remove these contributions. We will add the 
remaining subdominant discretization errors a posteriori 
to the systematic error budget after the chiral fit. 

In addition to the pion masses and pion/kaon ener¬ 
gies, several parameters enter the expressions in Eqs. (31) 
and (32). For completeness, we compile the values of the 
fixed parameters in our chiral fits in Table V. We use the 
lattice spacings and low-energy constant fi obtained in 
Ref. [16] from the RBC/UKQCD analysis of light pseu¬ 
doscalar meson masses and decay constants. We use 
the PDG value of = 130.4(2) MeV [10], and take 
A x = 1 GeV for the scale in the chiral logarithms. We 
use the B*Bn coupling constant gb = 0.57(8) obtained 
in our companion analysis also using the RBC/UKQCD 
domain-wall+Iwasaki ensembles and the RHQ action for 
the 5-quarks [51]. 

We perform correlated chiral-continuum fits to the 
data calculated on all five sea-quark ensembles listed in 
Table I using the full-QCD NLO SU(2) hard-pion/kaon 
HMyPT expressions. For B —$■ niv, we include dis¬ 
crete lattice momenta up to p"/ = 27r(l, 1,1)/L, which 


1 Recent three- and four-flavor lattice-QCD calculations typically 

give values for in the range of about 300-400 MeV [46-50]. 

The 2013 Flavor Lattice Averaging Group (FLAG) review quotes 
the range = 339(17) MeV for three active flavors [11]. To 
be conservative, we take a slightly larger value Aqcd = 500 MeV 
for the power-counting estimates throughout this work. 
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corresponds to ss 0.78 GeV on the coarser ensembles. 
For B s —» Klv, where the statistical errors are smaller, 
we include momenta up to pk = 27r(2, 0, 0 )/L, or ss 
0.91 GeV on the coarser ensembles. For the pion masses 
in Eqs. (31) and (32), we use the tree-level expression 
in Eq. (35). We obtain the physical form factors af¬ 
ter the chiral-continuum fit by setting the light quark 
mass to the physical average u/d- quark mass a^m u d = 
0.00102(5) [16] and the lattice spacing to zero. 

Figure 6 shows the resulting fits, which all have good 
X 2 /dof and p- values. We do not observe any statistically- 
significant lattice-spacing dependence for any of the form 
factors, and cannot resolve the coefficients of the a 2 terms 
in the four fits. Dropping the a 2 term altogether does not 
reduce the fit quality, and we consider this alternate fit as 
one of many possibilities when estimating the systematic 
uncertainty due to the chiral-continuum extrapolation in 
Sec. IV A. We observe a mild sea-quark mass dependence 
for f^ sK j and cannot resolve any sea-quark mass depen¬ 
dence in the other form factors. Dropping the term pro¬ 
portional to M 2 reduces the p- value of the f^’ K to ~ 5%, 
which is still acceptable, and does not impact the quality 
of the other fits. Again, we consider this alternative when 
estimating the chiral-continuum extrapolation error. Fi¬ 
nally, we do not see any evidence for the onset of chiral 
logarithms given that our lightest pion « 290 MeV 
is still quite heavy, and consider fits without the loga¬ 
rithms in Eqs. (31) and (32) among the alternate fits for 
assessing the systematic uncertainty. 

As a consistency check of our chiral-continuum extrap¬ 
olation, we can use our B — > 7 t£v form-factor fit results 
to obtain a rough estimate for the B*Bn coupling at 
lowest order in the l/m b expansion of HMyPT. From 
our preferred fits of / 1 | and f±, we find that the ratio of 
leading-order coefficients gives 

g b » cW/cW = 0.35(15) (36) 

where the error is statistical only (and does not include 
omitted higher-order corrections in the chiral and 1 /m b 
expansions). The value for g b in Eq. (36) is consistent 
with our independent determination of the B*Btt cou¬ 
pling in Ref. [51], and mostly independent of the input 
value of g b in the chiral logarithms. 

We also considered chiral-continuum extrapolations of 
the B —> 7r tv form factors using NLO SU(3) HMyPT, in 
which the logarithms have explicit strange-quark mass 
dependence, but were unable to obtain good fits for /^ 7r . 
Fits of to NLO SU(2) “soft-pion” yPT, in which the 
logarithms have explicit dependence on the pion energy, 
also failed to describe the data. All fits tried led to ac¬ 
ceptable p- values for the case of due to the fact that 

the shape is largely dictated by the B* pole term in the 
denominator. Finally, we tried supplementing the NLO 
expressions for the B —> i t£u form factors with NNLO 
analytic terms. The resulting partly NNLO fits yielded 
form-factor results consistent with those from our pre¬ 
ferred fits, but with significantly larger uncertainties due 


to the fact our data could not resolve any of the higher- 
order terms. 


IV. ESTIMATION OF SYSTEMATIC ERRORS 

We now discuss the sources of systematic uncertainty 
in our determinations of the B —> i\iv and B s —» Klv 
form factors. Each uncertainty is discussed in a sepa¬ 
rate subsection. We visually summarize the error bud¬ 
gets for the form factors versus q 2 in Fig. 7, and provide 
a detailed numerical error budget for the form factors at 
three representative q 2 values within the range of simu¬ 
lated lattice momenta in Table VI. The form factors at 
these three points will be used later in Sec. V for the 
extrapolation to q 2 = 0 via the z expansion. 

In cases where the estimation of a systematic uncer¬ 
tainty requires the explicit variation of simulation param¬ 
eters, we use the a « 0.11 fm ensemble with ami = 0.005, 
and take the dependence of that ensemble to be repre¬ 
sentative of all ensembles. We choose this ensemble be¬ 
cause it has very high statistics, and therefore allows us 
to most reliably measure the dependence of the form fac¬ 
tors on on the input parameters. We expect the behavior 
of the form factors on this ensemble to provide conserva¬ 
tive bounds on the errors since it has the largest lattice 
spacing and heaviest kaons. 

A. Chiral-continuum extrapolation 

We estimate the systematic uncertainty due to the 
chiral-continuum extrapolation of the B —> 7r and 
B s —► K form factors by varying the chiral-continuum 
fit Ansatze. We consider the following fit alternatives: 

• standard HMyPT including explicit Ep depen¬ 
dence in the chiral logarithms 

• omitting the term proportional to a 2 in Eqs. (31) 
and (32) 

• omitting the term proportional to M 2 in Eqs. (31) 
and (32) 

• omitting terms proportional to a 2 and M 2 in 
Eqs. (31) and (32) 

• analytic fits omitting the chiral logarithms in 
Eqs. (31) and (32) 

• analytic fits omitting the chiral logarithms and the 
term proportional to a 2 in Eqs. (31) and (32) 

• varying the value of in the coefficients of the 
chiral logarithms from /o= 112(2) MeV [16] in the 
chiral limit to Jk = 156.1(8) MeV [10] 

• varying the B* Bn coupling in the coefficients of 
the chiral logarithms g b = 0.57(8) by plus/minus 
one standard deviation [51] 
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FIG. 6. Chiral-continuum extrapolation of the B —» irlv (upper plots) and B s —> Klv (lower plots) form factors from correlated 
fits using NLO SU(2) hard-pion/kaon HM\PT. Fits of f± are on the left and of fu are on the right. In each plot, the colors 
distinguish between data points on the five different ensembles: circles and squares correspond to the a ss 0.11 fm data and 
triangles and diamonds to a ~ 0.086 fm data. The colored fit curves show the interpolation/extrapolation in pion/kaon energy: 
the fit function is evaluated at the unphysical sea-quark masses and nonzero lattice spacings on the different ensembles, such 
that the curves should go through the data points of the same color. The continuum, physical-quark-mass form factors are 
shown as a function of pion/kaon energy by the black lines with gray error band. The vertical dashed line on the left-hand side 
of each plot shows the physical pion or kaon mass. 


• varying the scalar pole mass Mb*( 0 + ) = 5.63 GeV 
in fQ sK by plus/minus 100 MeV 

• omitting the data point at zero momentum 

• omitting the data point at the highest momentum 
p = 2n/L(2, 0,0) for/®;* 

• excluding ensembles with pion masses > 400 
MeV 

Figure 8 shows the relative changes of the form-factor 
central values under each fit variation 

A/ i = i/r ef --/f lt i//r ef -. (37) 

where i = {0, +}. We take the largest difference between 
our preferred fit and any of the alternate fits as system¬ 
atic uncertainty due to the chiral-continuum extrapola¬ 
tion. We do not use fits with p-values below 1% or those 


that cannot resolve the coefficients within statistical un¬ 
certainties for our error estimate. Thus we exclude the 
fit omitting ensembles with pion masses M n > 400 MeV 
and the fit using standard “soft-pion” HMyPT. 

For each form factor, we obtain the largest difference 
from our preferred fit using the following variation: 

f^ : analytic for 18.7 GeV 2 < q 2 < 22.7 GeV 2 
and omitting the term elsewhere, 
fo* : analytic, 
f+ sK : analytic, 

: omitting the a 2 and M v terms, 

We therefore use these fits to obtain the q 2 -dependent 
chiral-continuum extrapolation errors quoted in Table VI 
and shown in Fig. 7. 
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TABLE VI. Error budgets for the B —» n£u and B s —> Klv form factors at three representative q 2 values in the range of 
simulated lattice momenta. For convenience, we also show the corresponding pion or kaon energy, Ep. Errors are given in %. 
The total error is obtained by adding the individual errors in quadrature. 



/r 

/<?" 

rB s K 

/+ 

rB s K 

JO 

E P [GeV] 

0.85 

0.50 

0.27 

0.85 

0.50 

0.27 

1.07 

0.77 

0.53 

1.07 

0.77 

0.53 

q 2 [GeV 2 ] 

19.0 

22.6 

25.1 

19.0 

22.6 

25.1 

17.6 

20.8 

23.4 

17.6 

20.8 

23.4 

/(«") 

1.21 

2.27 

4.11 

0.46 

0.68 

0.92 

0.99 

1.64 

2.77 

0.48 

0.63 

0.81 

Statistics 

7.9 

5.9 

12.4 

7.3 

4.6 

3.3 

4.1 

3.4 

3.2 

3.4 

2.7 

2.6 

Chiral-continuum extrapolation 

6.3 

5.0 

6.2 

10.9 

7.6 

5.8 

3.2 

2.8 

2.5 

5.0 

4.9 

5.1 

Light-quark mass m u d 

0.3 

0.2 

0.2 

0.4 

0.3 

0.2 

0.1 

0.1 

0.1 

0.1 

0.1 

0.1 

Strange-quark mass m s 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.1 

0.1 

0.1 

0.0 

0.0 

0.0 

Lattice-scale uncertainty 

2.0 

2.0 

2.0 

2.2 

2.2 

2.2 

2.0 

2.0 

2.0 

2.2 

2.2 

2.2 

RHQ parameter tuning 

0.9 

0.9 

0.8 

1.0 

1.0 

1.0 

0.9 

0.9 

0.9 

1.0 

1.0 

1.0 

Renormalization factor 

0.8 

0.8 

0.7 

1.6 

1.6 

1.7 

0.9 

0.8 

0.8 

1.6 

1.6 

1.7 

Finite volume 

0.5 

0.4 

0.3 

0.7 

0.5 

0.4 

0.2 

0.2 

0.2 

0.2 

0.1 

0.1 

Heavy-quark discretization errors 

1.8 

1.8 

1.8 

1.7 

1.7 

1.7 

1.8 

1.8 

1.8 

1.7 

1.7 

1.7 

Light-quark & gluon discretization errors 

1.1 

1.1 

1.1 

1.1 

1.1 

1.1 

1.3 

1.3 

1.3 

1.3 

1.3 

1.3 

Isospin breaking 

0.7 

0.7 

0.7 

0.7 

0.7 

0.7 

0.7 

0.7 

0.7 

0.7 

0.7 

0.7 

Total (%) 

10.6 

8.4 

14.3 

13.6 

9.6 

7.6 

6.2 

5.5 

5.3 

7.1 

6.7 

6.8 


B. Lattice-scale uncertainty 

We tuned the parameters of the 6-quark action to re¬ 
produce the experimental value of the B s meson mass, 
and carry out our analysis in terms of dimensionless ra¬ 
tios over ihg j to remove all explicit dependence on the 
lattice scale. We then obtain the form factors and mo¬ 
mentum transfers in physical units by multiplying by the 
appropriate power of Mb, = 5.367(4) GeV. The uncer¬ 
tainties in the form factors due to to the experimental 
error on Mb, are negligible. 

We do, however, still need to consider the implicit de¬ 
pendence on the lattice spacing through the parameters 
of the 6-quark action. We estimate the size of this depen¬ 
dence, by computing the form factors f± and f\\ for seven 
sets of RHQ parameters. We then calculate the slopes 
with respect to the parameters - A//Amoa, Af / Acp 
and A//A£ - for all momenta used in the analysis. Next 
we multiply each slope by the uncertainty in the corre¬ 
sponding RHQ parameter due to the lattice spacing from 
Table II, e.g. Af/Amo x a^ nQa . Finally, we add the indi¬ 
vidual contributions from the three RHQ parameters in 
quadrature to obtain the total systematic error due the 
lattice spacing. 

We examined the slopes with respect to the RHQ pa¬ 
rameters for both B — > n and B s — > K , and found 
them to be consistent. We therefore base our estimates 
for the systematic uncertainty due to the lattice spac¬ 
ing on the slopes obtained for the B s —> I\ form fac¬ 
tors because the smaller statistical errors in B s —> K 
enable the slopes to be resolved more precisely. Figure 9 
shows the slopes of the B s —> Ktv form factors with 
respect to the {moa,cp,(} on the a ~ 0.11 frn ensem¬ 
ble with ami = 0.005. For this slope estimate, we use 
the unimproved heavy-light vector current from Eq. (10). 
We find the largest slopes at p = 27t( 2,0,0) for f± and 


p = 27r(l, 1, 0) for /y. Following the procedure outlined 
above, we estimate lattice-spacing errors in f± and /y 
of 1.9 % and 2.2 %, respectively. In the continuum this 
corresponds to errors on f + (/ 0 ) of 2.0% (2.2%) which 
we take for both B s —> K and B —> tt. 


C. Light- and strange-quark mass uncertainties 

Here we estimate the error in the form factors due to 
the uncertainty in the light-quark mass and the mistim¬ 
ing of the strange sea quark. For clarity we discuss sepa¬ 
rately each place where the light- or strange-quark mass 
enters the analysis. 


1 . u/d-quark mass uncertainty 

We obtain the physical form factors f± and /y af¬ 
ter the chiral-continuum fit by evaluating Eqs. (31) and 
(32) at the physical average u/d -quark mass a 32 m^j) ys = 
0.00102(5). We estimate the error in the form factors due 
to the light-quark mass uncertainty by varying mj^ y& by 
plus/minus one sigma. For B —> n the central value shifts 
by 0.2 — 0.3% for /+”' and 0.2 — 0.4% for /(f 71 ', while for 
B s —► K both f+’ K and fQ aK change by 0.1%. 


2. Strange sea-quark mistuning 

Our preferred chiral-continuum fit employs SU{2 ) chi¬ 
ral perturbation theory, in which the strange quark mass 
is integrated out, so our fit function has no explicit depen¬ 
dence on m s . Further, at each lattice spacing, results for 
the form factors are only available at a single value of the 
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FIG. 7. VisLialization of the error budgets for the B — > tt£i/ (upper plots) and B a — > Kiv (lower plots) form factors. Error 
budgets for f± are on the left and of fu are on the right. The curves from bottom-to-top show the increase in the total 
percentage error as we add each individual source of error in quadrature. In each plot, the left y-axis label shows the squared 
error, while the right y-axis label shows the error in the form factor. For readability, we have combined all of the sources of 
uncertainty that we estimate to be below ~ 1% into a single entry labeled “other systematics.” The three vertical lines in each 
plot show the location of the synthetic data points used in the subsequent extrapolation to q 2 = 0. Detailed error budgets at 
these q 2 values are given in Table VI. 


strange sea-quark mass, so we cannot directly compute 
the strange sea-quark mass dependence of /y and f±. We 
therefore study the light sea-quark mass dependence and 
use it to bound the strange sea-quark mass dependence. 
We cannot resolve any light sea-quark mass dependence 
within statistical uncertainties, and expect the strange 
sea-quark mass dependence to be even smaller. Thus we 
take the error due to mistuning the strange sea-quark 
mass to be negligible. 


3. Valence strange-quark mass uncertainty 

The B s —> K form factors have explicit strange 
valence-quark mass dependence. The strange-quark 
masses employed in our simulations differ slightly from 
the physical, tuned values a 2 4 mj /lys = 0.0379(11) and 
a32m.P hvs = 0.0280(7) [16]. To study the valence strange- 


quark mass dependence, we calculated the B s —> K form 
factors on the a ss 0.11 fm, ami = 0.005, ensemble with 
two additional spectator-quark masses of a 2 ifh s = 0.033 
and 0.043. Figure 10 shows the valence-quark mass 
dependence of the B s —>• K form factors; we observe 
the largest slopes for /y at p = (0,0,0) and for f± at 
p = (1,0,0). Multiplication of these measured slopes 
by the discrepancy between the simulated and tuned 
strange-quark masses, A(m s ) = (m s — mf ys ) = 0.004, 
leads to estimates for the error due to mistuning the va¬ 
lence strange-quark mass of about 0.1% for f + and below 
this for /o (which we consider as negligible). 


D. RHQ parameter uncertainty 

We compute the semileptonic form factors using the 
nonperturbatively tuned RHQ parameters obtained in 
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FIG. 8. Relative change of the form-factor central value under the considered fit variations for B —► wlv (upper) and B a —> Klu 
(lower). In each plot, the shaded band shows the statistical uncertainty of the preferred fit. The three vertical lines show the 
location of the synthetic data points used in the subsequent extrapolation to q 2 = 0. 


Ref. [18] and given in Table II. The RHQ parameters 
have four significant sources of uncertainty: statistics, 
heavy-quark discretization errors, lattice scale, and the 
experimental inputs. We already discussed the uncer¬ 
tainty due to the lattice scale in Sec. IV B. We follow 
the same approach for propagating the uncertainty in 
the RHQ parameters due to heavy-quark discretization 
errors and experimental inputs. We multiply the esti¬ 
mated slopes of the form factors with respect to changes 
in mo, cp, and £ (shown in Fig. 9) by the uncertainties 
in the corresponding parameters due to heavy-quark dis¬ 
cretization errors and experimental inputs. Adding the 
contributions from the three RHQ parameters and the 
two uncertainty sources in quadrature, we obtain error 
estimates for of 0.8-0.9%, of 1.0%, f+ sK of 
0.9%, and fo’ K of 1.0%. 

We neglect the statistical uncertainties in the RHQ 
parameters in our final analysis, after checking that they 
have a negligible impact on the form factors. On the 


a ~ 0.11 fm, ami = 0.005 ensemble, we computed the 
form factors with seven sets of RHQ parameter values. 
We then used the approach detailed in Refs. [14, 18] to 
interpolate to the tuned RHQ parameters. This proce¬ 
dure automatically propagates the statistical errors in 
the RHQ parameters via the jackknife. Indeed we find 
that the statistical errors obtained from the two proce¬ 
dures are identical. Thus we do not need to perform the 
more complicated and computationally expensive proce¬ 
dure of interpolating to the tuned RHQ parameters in 
our analysis. 


E. Heavy-quark discretization errors 

The RHQ action gives rise to nontrivial lattice-spacing 
dependence in the form factors in the region mqa ~ 1. 
To estimate the size of the resulting discretization errors, 
we use the same power-counting approach as in our com- 
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heH ~p K = 2 jt(0,0,0)/L ^ ~p K = 2 jt( 1,1,0)/L ~p K = 2jt(2,0,0)/L 

p K = 2 jt(1,0,0)/L hsh f K = 2 jz(1,1,1)/L 



FIG. 9. RHQ parameter dependence of the B s —» K form factors f± (left) and f\\ (right) on the 24 3 ensembles with ami = 0.005 
using the unimproved heavy-light vector current in Eq. (10). The slopes are normalized using the form factors obtained at the 
central set of RHQ parameters. From left to right, the plots show the dependence on mod, cp, and Q The colored lines show 
the results of a linear fit to the three data points at each momentum. The black vertical lines indicate the tuned values of the 
RHQ parameters. The shaded vertical bands indicate the systematic errors in the RHQ parameters due to the lattice-scale 
uncertainty. For clarity, data points at equal RHQ parameter values are plotted with a slight horizontal offset. 


panion papers on bottomonium masses and splittings [18] 
and !3( s )-meson decay constants [14]. 

We tune the parameters of the operators in the 
dimension-5 RHQ action nonperturbatively, such that 
the leading heavy-quark discretization errors from the 
action are of 0(a 2 ). We use an (D(a s a)-improved vec¬ 


tor current and calculate the improvement coefficient to 
1-loop; therefore the leading heavy-quark discretization 
errors from the current are of 0(a 2 a, a 2 ). Because we 
use the same actions and simulation parameters as in 
our earlier calculation of the -meson leptonic decay 
constants [14], the numerical error estimates are almost 
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TABLE VII. Percentage errors from mismatches in the action and current for the bottom quark on the 24 3 and 32 3 ensembles. 
For this estimate, we calculate the mismatch functions for the nonperturbatively tuned parameters of the RHQ action from 
Table II. We estimate the size of operators using HQET power counting with Aqcd = 500 MeV and the coupling constant 
af s (l/a) = 1/3 on the 24 3 ensemble and 0.22 on the 32 3 ensembles. To obtain the total, we add the individual errors in 
quadrature, including each contribution the number of times that operator occurs. Contribution E is counted twice, and 3 is 
counted twice for /m and four times for f±. The definitions of operators “E,” “Xi,” “ X 2 ,” “Y,” and “3” and expressions for 
the mismatch functions are given in Appendix B of Ref. [14] . 



<U(a 3 ) error 
from action 

0(a 2 ) errors 
from current 

(D(aja) error 
from current 

Total (%) 


E 

V, X 2 Y 

3 

/ll u 

a as 0.11 fm 

0.55 

0.67 1.27 1.34 

1.48 

2.97 3.64 

a 0.086 fm 

0.42 

0.46 0.85 0.91 

0.55 

1.65 1.82 


identical in the two works. The same operators con¬ 
tribute in both cases, but enter a different number of 
times for the spatial and temporal vector currents. Ta¬ 
ble VII quotes the estimate of heavy-quark discretization 
errors from the five different operators in the action and 
current on the 24 3 and 32 3 ensembles, and we refer the 
reader to Sec. V. E. and Appendix B of Ref. [14] for 
details. We take the size of heavy-quark discretization 
errors in our calculation of the B —> and B s — > K£u 

semileptonic form factors to be the estimate on our finer 
a -1 = 2.281 GeV lattices, which is ~ 1.7% for the lattice 
form factors f\\ and ~ 1.8% for f±. These lead to errors 
in the continuum form factors /+ and fo of ~ 1.8% and 
~ 1.7%, respectively. 


F. Light-quark discretization errors 

The dominant discretization errors from the light- 
quark and gluon sectors are of 0 ((oAqcd) 2 ) from the 
action, and is about 5% using Aqcd = 500 MeV. We 
remove this 0 ((oAqcd) 2 ) error by including a term pro¬ 
portional to a 2 in the chiral-continuum extrapolation 
(see Eqs. (31) and (32)). Then the leading light-quark 
and gluon discretization errors in the heavy-light vector 
current are of 0(a s am q , (arh q ) 2 , a 2 a Aqcd, ( ap ) 2 ) where 
m q denotes the bare lattice mass. The first entry of 
0(a s arh q ) leads to estimated errors of ~ 0.1% in B —► n 
and ~ 0.6% in B s —> K form factors on the 32 3 ensem¬ 
bles. The second is negligible (< 0.1%) and the third is 
estimated to be ~ 1.1% in both £? —» n and B s —» K 
form factors. Adding these contributions in quadrature, 
we estimate the total uncertainty from light-quark and 
gluon discretization errors in the heavy-light current to 
be 1.1% in B — y tt and 1.3% in B s —► K. 

We do not observe any evidence of sizable momentum- 
dependent discretization errors in our data. Figure 2 
shows that the pion and kaon energies and amplitudes 
are consistent with continuum expectations, and smaller 
than power-counting estimates of 0((ap) 2 ). Thus we 
do not include a systematic error due to momentum- 
dependent discretization errors. 


G. Renormalization factor 


We renormalize the lattice form factors using a mostly 
nonperturbative approach in which we separate Zy into 
three components. We consider the uncertainties from 
these three multiplicative factors separately, and then 
add them in quadrature to obtain the total error on the 
form factors. 

For Zy , we use the nonperturbatively determined value 
of the axial-current renormalization factor Za in the chi¬ 
ral limit from Ref. [16]. We can neglect the statistical 
uncertainty in Za (which is only 0.02% on the finer en¬ 
sembles) and the difference between Zy and Za (which 
is about 0(am les ) ~ 7 x 10 -4 at a « 0.086 fm). For 
Zy, we use the nonperturbative determination from [14]. 
The statistical uncertainty in Zy on the finer ensemble 
is 0.15%. We conservatively estimate the perturbative 
truncation error in py to be the full size of the 1-loop 
correction at the finer a ~ 0.086 fm lattice spacing, which 
leads to 1.7% for py 0 and 0.6% for pv, • These are sig¬ 
nificantly larger than what we would estimate for two- 
loop contributions from naive power counting. Taking 
a s ~ 0.23 on the coarser lattice spacing and a coefficient 
of 1/(47 t) 2 from two loop-suppression factors, we would 
obtain an estimate of 0.03%. Even with a coefficient 
of 1 / 7r 2 , we would obtain an estimate of 0.5%, which is 
slightly smaller than the perturbative uncertainty that 
we assign to pv t and 3 times smaller than the error we 
assign to pv a ■ Because we use the values of pv fl and Zy in 
the chiral limit, we must consider the errors due to the 
nonzero physical up, down, and strange-quark masses. 
The leading quark-mass dependent errors in pv and Zy 
are 0{a. s afh q ) and 0((arh q ) 2 ), respectively, but these are 
already accounted for in our estimate of light-quark and 
gluon discretization errors (see Sec. IV F). Thus we do 
not count them again here. 

Perturbative truncation errors are by far the dominant 
source of uncertainty in the renormalization factor, and 
the quadrature sum of the three error contributions is 
1.7% for /y and 0.6% for f±. 
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k* f K = 2 jt( 0,0,0)/L -p K - 2 ji( 1,1,1 )/L 

heh J? K - 2ji(1,0,0)/L ~p K = 2jt(2,0,0)/L 

^ = 2jt( 1,1,0)/L 



a 24 ( m s + m res ) 



FIG. 10. Valence strange-quark mass dependence of the 
B s —» K form factors f± (top) and /n (bottom) on the 
a « 0.11 fm ensemble with ami = 0.005. The slopes are 
normalized by the form factors obtained with the strange- 
quark mass used in our production simulations. The colored 
lines show the results of a linear fit to the three data points 
at each momentum. The black vertical line with error band 
shows the total (statistical plus systematic) uncertainty in the 
physical strange-quark mass [16]. For clarity, data points at 
equal strange-quark masses are plotted with a slight horizon¬ 
tal offset. 


H. Finite-volume errors 


We compute the form factors on a finite-sized lattice. 
We estimate the effect of the finite spatial volume using 
one-loop finite-volume SU(2) hard-pion yPT, in which 
loop integrals are replaced by a sum over lattice sites. 
Only a single integral enters the NLO SU(2) hard-pion 
ChPT expression. The correction to Eqs. (33) and (34) 
to account for the finite spatial volume is given by a sum 


over modified Bessel functions [52, 53]: 



Ml log 


M 2 \ 4 ^ K^M^L) 

A 2 / + L ^ |f| 

7 r^O 1 1 


(38) 

From Eq. (38), the ensemble with the lightest quark mass 
receives the largest correction. For B —> tt we find correc¬ 
tions to /j_ (/y) of 0.3-0.4% (0.6-0.8%), while for B s —> I\ 
corrections to of 0.2-0.3% (0.4-0.5%). These re¬ 

sult in the following errors on the continuum form factors: 
0.3-0.5% for /fV 0.4-0.7% for /, 0.2% for and 

0.1-0.2% for f*' K . 


I. Isospin breaking 


Our B —> irlv and B s Klv form factors are cal¬ 
culated in the isospin limit. The form factors of the 
charged and neutral B (B s )~ mesons, however, differ due 
to both the masses and the charges of the constituent 
light u and d quarks. The leading quark-mass contri¬ 
bution to the isospin breaking from the valence-quark 
masses is of 0((md — m u )/AQcn) ~ 0.5%, which is 
obtained using the determination of the quark masses 
(wd — m u ) = 2.35(8)(24) MeV from FLAG [11] and 
Aqcd = 500 MeV. The difference between the it- and 
d-quark masses in the sea sector should have a negligible 
effect on the B —> Triv ( B s —► K£u) form factor because 
the sea quarks couple to the valence quarks through 1 = 0 
gluon exchange, and they give only the uncertainty of 
0(((rrid — m u )/ Aqcd) 2 ) ~ 0.003%. The electromagnetic 
contribution to the isospin breaking is expected to be 
0( a s) ~ 1/137 ~ 0.7% which is the typical size of 1-loop 
QED corrections. We therefore take 0.7% as the uncer¬ 
tainties due to the isospin breaking and electromagnetism 
effects. 


J. Correlation matrices 

In the next section we fit synthetic lattice data gener¬ 
ated at three values of q 2 to the ^-expansion to extend it 
to the full kinematic range. Thus, in addition to the sys¬ 
tematic uncertainties on the individual q ,2 -bins, we also 
need the correlations between q 2 values. Although it is 
straightforward to obtain the statistical correlations fur¬ 
ther explanation is needed for the systematic error cor¬ 
relations. 

The chiral-continuunr extrapolation error is estimated 
by varying the fit function and parametric inputs. This 
procedure does not provide any information on correla¬ 
tions of the resulting systematic error between different 
q 2 -bins. Alternate chiral-continuum fits to our data with 
different fit functions do, however, exhibit highly sim¬ 
ilar statistical correlations between < 7 2 -bins. Hence we 
take the (normalized) statistical correlation matrix from 
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TABLE VIII. Normalized statistical (upper) and systematic 
(lower) correlation matrices for the B -A ir£i/ form factors at 
three representative q 2 values. 


q 2 [GeV 2 ] 

19.0 

f+ n 

22.6 

25.1 

19.0 

fo n 

22.6 

25.1 


19.0 

1.000 

0.868 

0.045 

0.663 

0.586 

0.541 

/+* 

22.6 

0.868 

1.000 

0.239 

0.591 

0.654 

0.616 


25.1 

0.045 

0.239 

1.000 

0.176 

0.188 

0.283 


19.0 

0.663 

0.591 

0.176 

1.000 

0.822 

0.836 

fo w 

22.6 

0.586 

0.654 

0.188 

0.822 

1.000 

0.941 


25.1 

0.541 

0.616 

0.283 

0.836 

0.941 

1.000 





/+" 



fo * 



q 2 [GeV 2 ] 

19.0 

22.6 

25.1 

19.0 

22.6 

25.1 


19.0 

1.000 

0.897 

0.245 

0.702 

0.663 

0.645 

}+ n 

22.6 

0.897 

1.000 

0.427 

0.639 

0.725 

0.719 


25.1 

0.245 

0.427 

1.000 

0.289 

0.342 

0.448 


19.0 

0.702 

0.639 

0.289 

1.000 

0.840 

0.840 

fo W 

22.6 

0.663 

0.725 

0.342 

0.840 

1.000 

0.948 


25.1 

0.645 

0.719 

0.448 

0.840 

0.948 

1.000 


TABLE IX. Normalized statistical (upper) and systematic 
(lower) correlation matrices for the B a —> Kiv form factors 
at three representative q 2 values. 


q 2 [GeV 2 ] 

17.6 

rB s K 

J+ 

20.8 

23.4 

17.6 

rB s K 

Jo 

20.8 

23.4 


17.6 

1.000 

0.868 

0.828 

0.799 

0.754 

0.702 

rB s K 

J+ 

20.8 

0.868 

1.000 

0.783 

0.677 

0.799 

0.764 


23.4 

0.828 

0.783 

1.000 

0.615 

0.703 

0.708 


17.6 

0.799 

0.677 

0.615 

1.000 

0.828 

0.755 

rB s K 

Jo 

20.8 

0.754 

0.799 

0.703 

0.828 

1.000 

0.974 


23.4 

0.702 

0.764 

0.708 

0.755 

0.974 

1.000 





rB s K 

J+ 



rB s K 

Jo 



q 2 [GeV 2 ] 

17.6 

20.8 

23.4 

17.6 

20.8 

23.4 


17.6 

1.000 

0.939 

0.921 

0.865 

0.843 

0.808 

rB s K 

J+ 

20.8 

0.939 

1.000 

0.913 

0.794 

0.860 

0.835 


23.4 

0.921 

0.914 

1.000 

0.760 

0.806 

0.801 


17.6 

0.865 

0.794 

0.760 

1.000 

0.889 

0.840 

rB s K 

Jo 

20.8 

0.843 

0.860 

0.806 

0.889 

1.000 

0.983 


23.4 

0.808 

0.835 

0.801 

0.840 

0.983 

1.000 


our preferred fit and multiply it by the estimated chiral- 
continuum extrapolation error at each q 2 value (For off- 
diagonal elements of the correlation matrix we use the 
product a q 2 cr q 2 .) We follow the same procedure to esti¬ 
mate the correlations between the (^-dependent system¬ 
atic error due to the light-quark mass uncertainty. We 
take the remaining systematic errors for which we do not 
assume any q 2 dependence to be 100% correlated. 

Tables VIII and IX present the normalized statistical 
and systematic correlation matrices, which enable the full 
reconstruction of the total covariance matrices using the 
values for B —> tt£i/ and B s —► K£v form factors and 
their errors from Table VI. 


V. FORM-FACTOR RESULTS 

In this section we extrapolate our B —> tt ( B s —> K) 
form-factor results from large q 2 , where we have our 
(synthetic) data, to q 2 = 0 using a model-independent 
parametrization based on the general properties of an- 
alyticity, unitarity, and crossing symmetry. We first 
give the expressions for the z-parametrizations used in 
our analysis in Sec. V A; we use the parametrization of 
Bourrely, Caprini, and Lellouch (BCL) for our preferred 
results, but also consider the parametrization of Boyd, 
Grinstein, and Lebed (BGL) as a cross-check. Next, 
in Sec. VB we extrapolate our synthetic lattice data to 
q 2 = 0; we present our preferred results for f + and /o in 
Tables XI and XII as coefficients of the z-expansion and 
the matrix of correlations between them. 

Use of the z-parametrization to describe semileptonic 
form factors has several advantages over other functional 
forms used in the literature [54, 55]. Because the abso¬ 
lute value of |z| is small in the semileptonic region, and 
the z-coefficients are constrained to be small by unitarity 
and heavy-quark symmetry, one needs only the first few 
terms in the expansion to accurately describe the form 
factor shape with a negligible truncation error. More¬ 
over, as the precisions of both the lattice calculations 
and experimental measurements improve, one may easily 
include higher-order terms in z as needed. Finally, com¬ 
parisons of the z-expansion parameters resulting from 
fits to different theoretical or experimental data sets en¬ 
able a meaningful quantitative comparison of the shapes, 
while a combined fit to lattice and experimental data en¬ 
ables a clean determination of | |. The z-expansion 

has therefore been adopted as the preferred method for 
obtaining exclusive V u b by experimentalists on Babar 
and Belle, the Heavy Flavor Averaging Group, and the 
Particle Data Group [2-5, 10, 56]. 


A. z-expansions of semileptonic form factors 


The B -A 7 t£v and B s -A KHv form factors are ana¬ 
lytic functions of q 2 except at physical poles and branch 
cuts above the production threshold. Therefore, given a 
suitable change of variables, they can be expressed as a 
convergent power series (see, e.g., [24, 25, 57-60]). Uni¬ 
tarity and heavy-quark power counting bound the size 
of the series coefficients. In the literature, the new vari¬ 
able is called z, and the class of functions are called z- 
parametrizations. Two such parametrizations commonly 
used to extrapolate the B -a ~kIv form factor are by 
Boyd, Grinstein, and Lebed (BGL) [24] and Bourrely, 
Caprini, and Lellouch (BCL) [25]. 

The change of variables from q 2 to z is given by 


, 2 , \ = - q 2 /t+ - ^1 ~ t 0 /t+ 

[q,o) + 


(39) 


where + M T ) 2 and f_=(M S(s) — M^) 2 ■ This 
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transformation maps the semileptonic region 0 < q 2 < f _ 
onto a unit circle in the complex z plane. The B —> nlv 
and B s —> K£is form factors can then be expanded as a 
simple power series in z: 


BCL use this constraint on the derivative of the form fac¬ 
tor to remove an independent degree of freedom from the 
series expansion in z. Thus they arrive at the following 
parametrization for the vector form factor: 


Pi(q 2 )M(I 2 ,to)fi(q 2 ) ='^2 i af\to)z(q 2 ,t Q ) k , (40) 

k -0 

where i = { 0 ,+} for the scalar and vector form factors, 
respectively. The free parameter to in Eq. (39) deter¬ 
mines the range of \z\ in the semileptonic region, and 
hence can be chosen to accelerate the series convergence. 
The “Blaschke factor” Pi(q 2 ) must be chosen to vanish 
at any subthreshold poles to preserve the correct ana¬ 
lytic structure of fi(q 2 ). For /+, the relevant state is 
the J p = 1” meson, while for /o, the relevant state is 
the J p = 0 + meson. As discussed earlier in Sec. Ill C, 
the scalar B* meson has not been observed experimen¬ 
tally, but is predicted to have a mass well above the Bn 
production threshold [44, 45]. Thus the functions Pi for 
B — ► niv are typically taken in the literature to be 

P+* (q 2 ) = z(q 2 ,M B *), P^(q 2 ) = \. (41) 

Finally, the outer function (f>i(q 2 . to) can be any analytic 
function of q 2 ; different choices for <pi correspond to dif¬ 
ferent z-parametrizations. 

The form factors that describe B —► nlv ( B s — > Klv) 
in the range 0 < q 2 < f_, when analytically continued, 
also describe Bn ( BK ) production for q 2 > t + . The 
coefficients of the z-expansion are therefore bounded by 
the fact that the rate of production of Bn (BK) states 
is less than the production rate of all states coupling to 
the b — > u current. In Ref. [24], Boyd, Grinstein, and 
Lebed choose the outer function (f>i so that the unitarity 
constraint on the series coefficients takes a particularly 
simple form: 


N 2 

E (“1*0 £ !. (42) 


where this holds for any value of N. The explicit func¬ 
tions for </>® GL and <^q GL and their numerical values can 
be found in Ref. [60]. When using the BGL parametriza¬ 
tion for subsequent z-fits, we use to = 0.65 1- as in 
Ref. [60], such that —0.341 < z < 0.216 (—0.144 < z < 
0.148) for B —> nils (B s —► Klv) decay. 

In Ref. [25], Bourrely, Caprini, and Lellouch (who only 
discuss / + ) choose a simpler outer function </>g CL = 
1. They also point out that the BGL form-factor 
parametrization does not obey the known asymptotic be¬ 
havior near the Bn production threshold lmf + (q 2 ) ~ 

( q 2 — t+)' (which is due to angular momentum conser¬ 
vation). Therefore, at q 2 = t + (z = —1), the derivative 
of the form factor must satisfy 


' df+ 

dz 


= 0. 


U(<i 2 ) = 


1 — q 2 /r 


b B * 


K -1 

E 6 - 

k -0 


(k) 

+ 


z k _ ( _ 1)fc -X * * 

K 


(44) 


where we label the BCL series coefficients bk to distin¬ 
guish them from the BGL coefficients afc. There is no 
analogous constraint to Eq. (44) on the value or deriva¬ 
tive of /o at any z, so one cannot remove a further degree 
of freedom in the series expansion for the scalar form fac¬ 
tor. We therefore use the following functional forms for 
the scalar form factors: 


K -1 


/T(9 2 ) = E b ° )zk > 


k =0 


K -1 


fo BsK (<l 2 ) = 


1 - q 2 /m 2 B , 


E b o zk ’ 


k=0 


(45) 

(46) 


where we include a pole at the theoretically predicted 
value M b .( 0 + ) = 5.63 GeV for B s —► Klv [44]. Equa¬ 
tion (45) has been called the “simplified series expan¬ 
sion” in the literature [61]. To minimize the error from 
truncating the z-expansion for the B —> nils form factor, 
BCL choose t 0 = t opt = (M B + M„) (\JM B - y/M ^) 2 , 
such that the magnitude of |z| < 0.280 is minimized in 
the semileptonic range. With the analogous choice for 
B s —► K£is, \z\ < 0.146 for the semileptonic range. 

Although the functional form of the BCL parametriza¬ 
tion is simpler than that of BGL, the unitarity constraint 
on the coefficients is more complicated [25]: 


E B 3k(to)by\to)b?\to)<l, (47) 

j,k=o 

OO 

Pjk (to) ^ ^ Vn(d^)dn+\j — (^o)? (48) 

n =0 


where rji is the Taylor coefficients in the expansion of the 
outer function 


Mb 


*(*) = , Mq 2 (z),t 0 ) {1 ; )2(1 , (49) 

4 (t + -t 0 ) (1 -zz*R 


z* = z(M|,t 0 ), 


(50) 


around z = 0. The values of Bjk for the B —> n£is and 
B s —» K£is form factors with the choice to = f op t are 
given in Table X. 

For the B —> n£is vector form factor, Becher and 
Hill [62] use heavy-quark power counting to provide an 
estimate for the sum of the coefficients: 



k—0 



z =-1 


(43) 


(51) 
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TABLE X. Matrix elements Bjk(to) that enter the unitarity 
bound on the BCL series coefficients for the choice to = t op t. 
The remaining coefficients can be obtained from the relations 
Bjti+k) = Bok and the symmetry property Bjk = Bkj. To 
derive these results we use the outer functions <f> + and <j >o in 
Eq. (7) of Ref. [60] with inputs from Ref. [25], giving x+ = 
5.03 x 10“ 4 and Xq 0) = 1-46 x 10~ 2 . 



Boo 

-Boi 

Bq2 

B 03 

Bo4 

Bo5 

/+" 

0.0197 

0.0042 

-0.0109 

-0.0059 

-0.0002 

0.0012 

fo Bw 

0.1062 

0.0420 

-0.0368 

-0.0406 

-0.0201 

-0.0057 

rti s K. 
r^sK 

Jo 

0.0115 

0.0004 

-0.0076 

-0.0007 

0.0018 

0.0004 

0.0926 

0.0137 

-0.0484 

-0.0174 

-0.0003 

0.0024 


where A is a typical hadronic scale. Taking A ~ 
1000 MeV, this would imply )G a 2 ~ 0.01, which is well 
below the bound from unitarity. Experimental measure¬ 
ments [2-5] and previous lattice calculations [27] con¬ 
firm this expectation. This argument also applies to the 
B s —> Ktv vector form factor, where we emphasize that 
Eq. (51) is only a rough constraint due to the imprecise 
scale A and omitted higher-order corrections in the OPE 
and 1 /m;,. 

B. Extrapolation of lattice form factors to q 2 = 0 

We now extrapolate our results for the B —> ir£v 
and B s —>• Kh> form factors to q 2 = 0 using the z- 
expansion. We first generate synthetic data points in the 
range of simulated data from the output of the chiral- 
continuum extrapolation. Recall that the continuum, 
physical quark-mass form factors are obtained from fits 
to Eqs. (31) and (32) by fixing M 2 to the physical value 
and a 2 —> 0. After these replacements, the physical form 
factors depend upon three independent functions of the 
pion or kaon energy Ep. We therefore generate three syn¬ 
thetic data points each for fo and /+ in order to ensure 
that the covariance matrix is not singular. In anticipa¬ 
tion of the z-fit, we choose the points to be evenly spaced 
in z (rather than q 2 ). The q 2 values and error budgets 
for the synthetic lattice data are given in Table VI. 

We fit our synthetic lattice data for the B —>• 7 r Iv and 
B s —> Klv form factors including statistical and system¬ 
atic correlations between q 2 values. For our preferred 
fit we use the BCL parametrization with the kinematic 
constraint /+(0) = /o(0) and use the theoretical estimate 
from heavy-quark power counting to constrain the sum of 
the coefficients of the vector form factor via Bayesian pri¬ 
ors. We study the central values and errors of the series 
coefficients as a function of the truncation K such that 
our final form-factor results include the truncation error. 
The complete z-fit results are given in Appendix B. We 
also compare to results using the BGL parametrization 
as a check. 

We first perform separate fits of f + and f 0 without 
imposing any constraints on the sum of coefficients. The 


results for B —> tt£v are given in the top two panels of 
Table XIX, and for B s —> K£v in the upper two panels 
of Table XX. The separate fits of /+ and f 0 for I\ = 2, 3 
are shown in the left-hand plots of Fig. 11 for B —> niv 
(upper) and B s —► K£v (lower). The synthetic lattice 
data points are correlated, and one must include a term 
quadratic in z to obtain a good fit (recall that for f + the 
expression with K = 2 includes a term proportional to z 2 
that is related to the z° and z 1 terms). The normaliza¬ 
tions are well determined by the lattice data, with 
central values that are stable within errors when going 
from K = 2 to K = 3. This is important because the 
normalization of the vector form factor plays a key role 
in the determination of V u b (see Sec. VIA). We cannot 
go beyond K = 3 because we have only three synthetic 
data points. 

In the separate fits to /+ and fo with K = 3, the kine¬ 
matic constraint /+( 0 ) = /o( 0 ) is automatically satisfied 
within uncertainties, but with large errors. We can there¬ 
fore impose the kinematic constraint /+(0) = /o(0). The 
results of the combined fits are given in the third panels 
of Tables XIX and XX. As expected, the constrained fits 
with K = 2 for both f + and f 0 have poor p-values, but 
the remaining fits tried are all of good quality. Adding 
the kinematic constraint (and only considering the good 
fits) has little impact on the results for the normalizations 
and even on the slopes (fc^/fe- 0 ^). It reduces the error 

on the curvatures (frfV^) as compared to the separate 
fits with K = 3, however, and consequently improves the 
determination of f+{q 2 = 0 ). 

Even with the kinematic constraint, however, the 
slopes and curvatures of the form factors are still not well 
determined by the lattice data, with errors ranging from 
25% to as much as 300%. For all fits considered, the sum 
of the coefficients B jkbjb k satisfy the unitarity con¬ 
straint. Further, for / + , the sum ^ Bj k bjb k is also con¬ 
sistent with expectations from heavy-quark power count¬ 
ing, Eq. (51), but with large uncertainties. We can there¬ 
fore use theoretical guidance from heavy-quark power 
counting to further improve our lattice form-factor de¬ 
termination. Keeping the kinematic constraint, we also 
constrain the sum of the coefficients of the B —> i t£v 
and B s —> Klv vector form factors with Bayesian priors 
based on their estimated size from heavy-quark power 
counting. For the hadronic scale in the heavy-quark es¬ 
timate we take 1000 MeV, with a generous uncertainty 
of ±500 MeV. Thus for the prior central value we use 
B = 0.01, and for the Gaussian prior width we use 
<jb = 0.03. We implement the Bayesian fit by minimizing 
the augmented x 2 [63], 

Xaug = X 2 + Xprior ) (52) 

where 

Xprior = ( B ~Y1 B jkbjb k ) 2 /(J 2 B . (53) 

The results for different truncations K are given in the 
bottom panels of Tables XIX and XX. The inclusion of 
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FIG. 11. Fits of the B — > tt£i/ (upper plots) and B s —> Kiu (lower plots) lattice form factors to the ^-expansion versus 
truncation K without constraints (left) and with the kinematic and heavy-quark constraints (right). The black open symbols 
show the synthetic data points with statistical (inner) and statistical ® systematic (outer) error bars. The curves with colored 
bands show the fit results with errors for different truncations K. We do not show unconstrained fits with K = 4 in the 
left-hand plot because we only have three synthetic data points; the inclusion of the kinematic and heavy-quark constraints 
allows us to perform the K = 4 fit shown in the right-hand plot. In the right-hand plot, we do not show the I\ = 2 combined 
fit because of the poor fit quality. 


the heavy-quark constraint improves the determinations 
of the slopes and curvatures, and leads to a reduction in 
the absolute error on f? n (0) by about a factor of 2 for 
B —> n£u for K = 3. The improvement in the error on 
f+ ,K { 0) is smaller but non-negligible, about 25%. 

After implementing the heavy-quark constraint, we are 
able to include an additional parameter in our fits and 
can consider expansions with K = 4. This enables us 
to study the stability of the central values and errors of 
the parameters with truncation K, and thus assess the 
systematic uncertainty associated with truncating the 21- 
expansion. The central values and errors for the normal¬ 
izations and slopes are stable when increasing the trunca¬ 
tion from K = 3 to K = 4, in most cases changing only in 
the last decimal place (except for the slope of fa sK {q 2 ), 
for which the results are still consistent within uncer¬ 
tainties). The combined fits of /+ and fo imposing the 
kinematic and heavy-quark constraints are shown versus 
the truncation I\ in the right-hand plots of Fig. 11 for 


B —> •niv (upper) and B s —> Ktv (lower). The central 
fit curves for K = 3 and K = 4 lie almost on top of 
each other, while the widths of the error bands and the 
uncertainties in /+(0) increase only slightly in going to 
K = 4. Thus we conclude that the K = 3 constrained 
fit includes the systematic uncertainty due to truncating 
the series in z. 

We therefore take as our preferred fits for B —> n£v 
and B s —> Kiv the results from the fit with K = 3 
for both / + and /o including the kinematic and heavy- 
quark constraints. This is the highest truncation K for 
which we still have more data points than fit parameters, 
and the uncertainties are comparable to the K = 4 fits. 
Figure 12 shows our preferred fits for B —>• tt£v (upper 
plots) and B s —> Ktv (lower plots) plotted versus z (left) 
and versus q 2 right. 

As a cross-check, we compare our preferred fit using 
the BCL parametrization to the analogous fit (also im¬ 
posing the kinematic and heavy-quark constraints, and 
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FIG. 12. Preferred K = 3 fit of the B —> nlv (upper plots) and B s —> Klv (lower plots) lattice form factors to the ^-expansion 
including the kinematic and heavy-quark constraints versus z (left) and versus q 2 (right). The black open symbols show the 
synthetic data points with statistical (inner) and statistical ® systematic (outer) error bars. The solid curves with error bands 
show the fit results for f+(q 2 ) and fo(q 2 )- 


to the same order in z ) using the BGL parametrization. 
Figure 13 overlays the results of the BCL and BGL fits 
for B —> 7r tv (left) and B s —> Klv (right). The fits to the 
different series expansions are consistent, indicating that 
our quoted form-factor uncertainties encompass the error 
due to truncating the z-expansion. The error bands from 
the BCL fits are narrower because the BCL form for /+ 
relates the coefficient of highest-order term in z to the 
coefficients of the lower-order terms. 

Tables XI and XII present our final results for the 
B —> nlv and B s —> Klv form factors as coefficients of 
the z-expansion and the matrix of correlations between 
them. These results are model independent and valid 
over the entire semileptonic region of q 2 . As we illustrate 
in the next section, they can be used in combined fits 
with experimental data to obtain the CKM matrix cle¬ 
ment \V u b\, or to make predictions for Standard-Model 
observables for these decay processes. 

It is interesting to compare ratios of these form factors 
to predictions from approximate symmetries of QCD. In 
the SU(3) limit ( rrid = m s ), the form factors for B —> nlv 


and B s —> Klv should be identical. Thus the ratios 
RM 2 ) = /f*(<? 2 )//f% 2 ) - 1, for * = {+,0}, pro¬ 
vide a measure of SU(3)-breaking in B —> light semilep¬ 
tonic form factors. Figure 14, left, plots these ratios 
for the full kinematic range. The results for f + and 
/o are similar. The deviations from unity are consis¬ 
tent with expectations from simple power counting of 
(m B — m<j)/ Aqcd ~ 20 %, but with large uncertainties. 

At large recoil (low q 2 ) and in the heavy-quark symme¬ 
try limit (mb/ Aqcd —t oo), the B —>• nils and B s —>• 
processes are each described by a single independent form 
factor as follows [64]: 

f o( q2 )= mB ^~ q /+(g 2 )- (54) 

This expression reduces to the kinematic constraint 
f+(q 2 = 0) = f 0 (q 2 = 0) at q 2 = 0. Figure 14, right, 
plots the ratio fo(q 2 )/f+(q 2 ) for the full kinematic range. 
The results are similar for B —> iriv and B s —> Klv. 
They agree exactly with the prediction from Eq. (54) at 
q 2 = 0 by construction because we imposed the kinematic 
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FIG. 13. Comparison of fits to the B —» ir£i/ (left) and B a —> K£v (right) lattice form-factor data using the BCL and BGL 
parametrizations. The black open symbols show the synthetic data points with statistical (inner) and statistical © systematic 
(outer) error bars. The gray and red curves with error bands show the BCL and BGL fits, respectively. 


TABLE XI. Central values, errors, and correlation matrix 
for the parameters of our preferred fit to the B —> i\iv form 
factors f + and fo to the BCL z-parametrization in Eqs. (44) 
and (45) using t 0 = t opt = (Mb + M n ) - \fMpf in 

Eq. (39) and pole mass M B * = 5.3252(4) GeV in f+(q 2 ). 



value 

bf 

Correlation matrix 

fe (l) & (2) A0) Ai) 

°+ u + u 0 °0 

b {2) 

°0 


0.412(39) 

1.000 

0.337 

-0.076 

0.679 

0.045 

0.100 


-0.511(184) 

0.337 

1.000 

0.150 

0.222 

0.698 

0.581 


-0.524(612) 

-0.076 

0.150 

1.000 

0.029 

0.436 

0.659 

6 (0) 

°o 

0.520(60) 

0.679 

0.222 

0.029 

1.000 

-0.258 

-0.224 

b (1) 

°0 

-1.657(182) 

0.045 

0.698 

0.436 

-0.258 

1.000 

0.564 

°0 

2.146(682) 

0.100 

0.581 

0.659 

-0.224 

0.564 

1.000 


constraint in our preferred z-fit, but are consistent with 
heavy-quark-symmetry expectations throughout the low 
q 2 region. 


VI. PHENOMENOLOGICAL APPLICATIONS 

In this section we present two phenomenological appli¬ 
cations of our form-factor results. 

First, in Sec. VIA, we use our results for the B —► tt£v 
form factors to determine the CKM matrix element \V u b\. 
We fit recent experimental measurements of the B —> tt£v 
differential branching fraction to the z-parametrization 
to obtain the slope b^P /b^P and curvature b^P /bP\ Con¬ 
firming that the lattice and experimental shapes are in¬ 
deed consistent, we then perform a combined z-fit of our 
numerical B —> 7 x£v form-factor data with the experi¬ 
mental measurements to obtain a model-independent de¬ 
termination of \V u b\. This method can also be applied to 
the decay B s —>■ K£u , once it has been observed experi- 


TABLE XII. Central values, errors, and correlation matrix 
for the parameters of our preferred fit to the B s —> K£v form 
factors /+ and fo to the BCL z-parametrization in Eqs. (44) 

and (46) using t 0 = t opt = ( M Bs + M K ) ( \JM Ba - s/Mk^J 
in Eq. (39) and pole masses M B * = 5.3252(4) GeV and 
Mg*(0 + ) = 5.63 GeV in f+(q 2 ) and fo(q 2 ), respectively. 




Correlation matrix 


value 

bf 


& V 2) 

b (0) 

°0 

b^ 

°0 

l (2) 

°0 

6 V 0) 0.338(24) 

1.000 

0.255 

0.146 

0.873 

0.603 

0.423 

lyp -1.161(192) 

0.255 

1.000 

0.823 

0.311 

0.954 

0.770 

b (2) -0.458(1.009) 

0.146 

0.823 

1.000 

0.346 

1.060 

0.901 

5q 0) 0.210(17) 

0.873 

0.311 

0.346 

1.000 

0.556 

0.479 

b^ -0.169(202) 

0.603 

0.954 

1.060 

0.556 

1.000 

0.965 

b ( 0 2) -1.235(880) 

0.423 

0.770 

0.901 

0.479 

0.965 

1.000 


mentally, to provide an alternate determination of \V u b\- 

Next, in Sec. VIB, we make predictions for Standard- 
Model observables for the decay processes B —> iriis and 
B s — ► K£u for both £ = h,t final-state charged lep¬ 
tons. (Here we use /r to indicate both muon and electron 
final states, for which the Standard-Model predictions 
are indistinguishable at the current level of precision.) 
We show results for the differential branching fractions, 
forward-backward asymmetries, and fi/r ratios (which 
are independent of |I4f,|). We only calculate observables 
that depend upon |I4b| for B s —► I\£v decays, using the 
value determined previously in Sec. VIA. Once the ex¬ 
perimental error on the branching fraction is commen¬ 
surate with the theoretical form-factor uncertainties, our 
B s —>• K£v form-factor results will enable a sufficiently 
precise determination of V u t, to illuminate the discrep¬ 
ancy between \V u b\ from inclusive B —> X u £v and exclu¬ 
sive B — » tt£v semileptonic decays. 
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FIG. 14. Left: SU(3)-breaking ratios R+(q 2 ) (red) and Ro{q 2 ) (blue). The dashed horizontal lines show the expected size of 
SU(3) breaking from simple power counting. Right: ratios of scalar to vector form factor fo(q 2 )/f+(q 2 ) for B —¥ itlv (red) 
and B s —> Kiv (blue). The black diagonal line shows the prediction in the heavy-quark symmetry limit. Errors are shown in 
shaded bands. In the left plot, errors shown are statistical only. 


A. Determination of \V u b\ from B — > n£v 

For the determination of \V u b\, we include the two 
most recent experimental measurements from BaBar, 
which are the untagged 6-bin (“BaBar 2010”) and 12-bin 
(“BaBar 2012”) analyses in Refs. [2, 4], Because the 12- 
bin analysis uses more data and different candidate selec¬ 
tions and cuts than the 6-bin analysis, the statistical cor¬ 
relations between the two data sets are quite small, and 
we treat the two data sets as statistically uncorrelated. 
There are some correlations between the systematic un¬ 
certainties in the two analyses, but these are estimated 
to be sufficiently small that they have a tiny impact on 
V u b [65]. We therefore treat the two BaBar analyses as 
fully independent. We also include the two most recent 
experimental measurements from Belle, which are the un¬ 
tagged analysis in Ref. [3] (“Belle 2010”) and the full- 
reconstruction tagged analysis in Ref. [5] (“Belle 2013”). 
The tagged and untagged data sets have little overlap. 
Further, the dominant systematic error in the tagged 
analysis is from the uncertainty in the tagging calibra¬ 
tion, which is not present for the untagged analysis. Thus 
we treat the Belle tagged and untagged analyses as inde¬ 
pendent. The BaBar and Belle data sets are statistically 
independent. The only commonality to the BaBar and 
Belle analyses is the use of the same event generation [65] . 
Because the event generation is not a significant source of 
uncertainty in the analyses, we treat the systematic un¬ 
certainties as uncorrelated between the BaBar and Belle 
data sets. 

We first fit the experimental measurements to the 
BCL z-parametrization to obtain the shape parameters 
(6+V^+’ ) )' F° r these fits, we do not impose any constraint 
on the sum of the coefficients '^2B mn b rn b n . Fits with 


truncation K = 2 are sufficient to obtain good y 2 /dof 
for three of the four data sets, but we perform fits with 
K = 3 in order to enable comparison of both the slopes 
and curvatures with those of the form factor ob¬ 

tained in the previous section. The numerical results for 
the K = 3 fits to the individual experimental data sets, 
as well from a combined fit to all experiments, are given 
in Table XXI. The fit to the BaBar 2010 data set has 
a somewhat large x 2 /dof that stems from the highest 
q 2 bin, for which the error on the measured differential 
branching fraction is small but the central value is low 
with respect to the other points. The inconsistency of 
the BaBar 2010 data leads the fit to all four experimental 
measurements to have a somewhat low, but still reason¬ 
able, p-value of 5%. Figure 15 shows the constraints on 
the slope (b+V&+^) versus curvature (b^/b^) from the 
different experimental measurements, as well as from the 
combined fit to all four measurements. The three most 
recent measurements agree at the 2er level, but display 
some tension with the BaBar 2010 result. Combining the 
information from all four experimental analyses improves 
the determination of the shape parameters significantly. 

Because we do not impose any constraint on the sum 
of the coefficients ^ we can check to see 

whether the experimental data is compatible with ex¬ 
pectations from heavy-quark power counting for the size 
of the series coefficients. Taking the determination of 
|V u b| = 3.63(12) x 10 -3 from CKM unitarity [7], we find 
a value for ^ B mn b m b n ~ 0.02 from the fit to all experi¬ 
mental data. This is consistent with the prediction from 
Eq. (51) taking a reasonable value for the heavy-quark 
scale A ~ 1.1 GeV, and validates the prior central value 
and width that we used to constrain ]V] 7? mra &m^n in our 
preferred z-fit of the lattice form factors in the previous 
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FIG. 15. Shape parameters b^/b® and b^/b^ from A' = 3 
BCL fits to the B — > tv£v form factor f+ w {q 2 ) (filled ellipse) 
and from experimental measurements of the B —> n£v branch¬ 
ing fraction [2-5] (patterned and empty ellipses). The colored 
ellipses show the constraints from the individual experiments, 
while the black ellipse shows the constraint from all experi¬ 
ments. For each determination, the inner and outer contours 
show the 68% and 95% allowed conhdence limits, respectively. 


section. 

Finally, before we fit the experimental and lattice data 
together to obtain \V u b\, it is important to check that 
their shapes are consistent. Figure 15 also shows the 
determination of the slope and curvature from our calcu¬ 
lation of (q 2 ) (see Table XI). The shapes of the lattice 
form factors and the experimental data are in good agree¬ 
ment, but the shape (as well as the overall normalization) 
is determined more precisely by experiment. This sug¬ 
gests that the error on \V u b\ can be minimized by per¬ 
forming a combined fit to the lattice and experimental 
data, as we now show. 

Table XXII shows the results for the BCL coefficients 
and \V u b\ obtained from a combined fit of the experimen¬ 
tal measurements for the B —> 7r iv differential branching 
fraction and the lattice determination of the form fac¬ 
tor f^(q 2 ), leaving the relative normalization \V u b\ as a 
free parameter to be determined in the fit. As in the 
experiment-only 2 -fits above, we do not constrain the 
sum of the coefficients B m nbmb n - We present results 
from separate fits to each experimental data set, as well 
as from a fit including all experimental data. The results 
for |I4t| from fits to the different experimental data sets 
agree within about lcr, and the p-value of the I\ = 3 fit 
to all data is 6%. We also show results for truncations 
K = 3,4, 5 to study the uncertainty due to truncating 


the expansion in z. The errors on \V u b\ remain the same 
size as the number of fit parameters increase, and the 
central value for the fit including all experimental data 
is unchanged. We take our final result 

\V ub \ = 3.61(32) x 1(T 3 (55) 

from the fit to all experimental data with K = 3. The 
quoted error on V u i, is the total uncertainty, and in¬ 
cludes both the theoretical error from the form factor 
and the experimental error (as well as the uncertainty 
from truncating the 2 -expansion). Figure 16 shows the 
preferred K = 3 BCL 2 -fit used to obtain \V u b\ plotted 
as (1 — q 2 /m 2 B ,)f + (q 2 ) vs. 2 (left) and as AB/Aq 2 vs. 
q 2 (right). 

Although we cannot precisely disentangle the error 
contributions, we can estimate the contribution to the 
error on \V u b\ from the lattice form-factor determination. 
Our most precise synthetic data point has a total statis¬ 
tical plus systematic uncertainty of 8.4%. If we assume 
that this is the lattice contribution to the 8.9% error | V u b\ 
in Eq. (55), this suggests that the experimental error con¬ 
tribution is approximately 2.8%. 

The combined 2 -fit optimally combines the available 
information from lattice and experiment in a model- 
independent manner, thereby providing a determination 
of 114 i> I that is both reliable and precise. We can quan¬ 
tify this statement by comparing the error on |I4b| from 
the simultaneous fit to the error obtained from the pre¬ 
viously standard approach. One can determine |14ft| by 
relating the measured partial branching fraction in an 
interval [<?min>9max] to the normalized partial decay rate 
calculated from the form factor as follows: 

A ^(9min I 9max)/|Kx6| 2 = A)A^min, ?Lx), (56) 

where 

^2 

A C(<7min>9max) E ^ | dq 2 p%\f + (q 2 )\ 2 . (57) 

9min 

The momentum range q 2 > 16 GeV 2 is typically used for 
comparison with lattice-QCD calculations because it is 
directly accessible in simulations and avoids (or at least 
minimizes) the need for an extrapolation. From our pre¬ 
ferred BCL parametrization (c/. Table XI) we obtain 

AC Bt (16 GeV 2 , q 2 max ) = 1.77(34) ps" 1 . (58) 

Table XIII shows the determinations of |I4ft| obtained 
from combining A£(16 GeV 2 , q 2 mx ) above with the differ¬ 
ent experimental measurements of AB(16 GeV 2 ,<?4ax)j 
and with their average, via Eq. (56). The results agree 
with those from the simultaneous 2 -fits, but with larger 
errors. In particular, the error on \V u b\ obtained using 
the average A£>(16 GeV 2 ,g4ax) from all experiments is 
10.0%, to be compared with the 8.5% error from our 
combined 2 -fit. Separate 2 -fits to the lattice and exper¬ 
imental data lead to a similar error on \V u b\- The error 
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FIG. 16. Model-independent determination of \V u i,\ from a combined fit of experimental measurements of the B —> tv£i/ 
branching fraction [2-5] and our lattice result for the B —> tv tv form factor /+(g 2 ) to the BCL 2 parametrization, Eqs. (44) 
and (45), with K = 3. The left plot shows (1 — q 2 /m|.)/+(ij 2 ) vs. z (where the experimental data have been rescaled by the 
value of \V u b\ determined in the fit), while the right plot shows AB/Aq 2 vs. q 2 (where the lattice points have been rescaled by 
|K&|). In both plots, the filled black circles show the lattice data, while the open colored symbols show the experimental data. 
The black curve with gray error band shows the fit result. 


TABLE XIII. Determinations of \V u b\ from a comparison 
of the measured B —> ntv partial branching fractions with 
the normalized partial decay rate A£ s,r (16 GeV 2 , g/Mx) = 
1.77(34) calculated from our preferred BCL paramterization 
of the vector form factor f+^(q 2 )- 



A 13(16 GeV 2 ,g^ax) x 10 7 

|I4b| x 10 3 

All 

368(19) 

3.69(37) 

BaBar 2010 [2] 

319(34) 

3.44(38) 

BaBar 2012 [4] 

369(32) 

3.70(39) 

Belle 2010 [3] 

398(30) 

3.84(40) 

Belle 2013 [5] 

386(51) 

3.78(44) 


on the normalization of the form factor &o in Table XI is 
9.4%, while the error on the normalization of the experi¬ 
mental branching fraction from the K = 3 fit to all exper¬ 
imental data &o|Tub| is 2.2%. Adding these in quadrature 
leads to a total error of 9.7%. Thus we conclude that 
the combined 2 -fit of all lattice and experimental data is 
indeed the best approach for minimizing the uncertainty 
on \V ub \. 


B. Standard-Model predictions for B — ¥ tt£v and 
B s —> Ktv observables 

The Standard-Model differential decay rate for Br a ) —> 
Ptv is given in Eq. (1). Using the experimentally mea¬ 
sured lepton and meson masses [10], we obtain predic¬ 
tions for the differential decay rate divided by |I4b| 2 - 
These are plotted for the muon and r-lepton final states 


in Fig. 17, where we use “muon” to denote decays to ei¬ 
ther of the light charged leptons (i = /x, e) throughout 
this section. Integrating the differential decay rates over 
the kinematically-allowed q 2 range gives 2 

T(B -4 n^u)/\V ub \ 2 = 6.2(2.5) ps" 1 , (59) 

T(B -»■ TTTv)/\V ub \ 2 = 4.3(1.2) ps" 1 , (60) 

r(E s -4 Kfj,v)/\V ub \ 2 = 4.55(1.08) ps" 1 , (61) 

T(B S -»■ Ktv)/\V ub \ 2 = 3.52(0.60) ps" 1 , (62) 

with errors of about 25-40% and 15-30% for the /x and r 
final states, respectively. We also use the determination 
of \V ub \ from our calculation of the B —> tv tv form factors 
(Eq. (55)) to make predictions for the B s —> Ktv differ¬ 
ential branching fractions for £ = /i, r. These are plotted 
in Fig. 18. For comparison, we also show the prediction 
for dB/dq 2 using the determination of \V u b\ from inclu¬ 
sive B —► X u tv decay [66] . The form-factor uncertainties 
are sufficiently small for q 2 > 13 GeV 2 that, given an 
experimental measurement of the branching fraction in 
this region with commensurate precision, one can distin¬ 
guish between the curves corresponding to |Uub| eX ci. and 
I Tub | inci. ■ Thus we anticipate that B s —> Ktv semilep- 
tonic decay will eventually play an important role in ad¬ 
dressing the current “|14h| puzzle.” 

Semileptonic decays to r leptons may be particularly 
sensitive to new physics associated with electroweak sym- 


2 In practice, the full kinematic range may not be accessible ex¬ 

perimentally, in which case the limits of integration here and 

throughout this section will need to be changed accordingly. 
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FIG. 17. Standard-Model predictions for the differential decay rate divided by | V^.61 2 for —> Plv decays to muon (left) and 

r-lepton (right) final states using our form-factor determinations in Tables XI and XII. 



FIG. 18. Standard-Model predictions for the differential branching fraction for B a — > K/iu (left) and B s — > Ktv (right) 
using our determinations of the B a — > Klv form factors in Table XII. Each plot shows predictions for dB/dq 2 using our 
determination of \V u b\ in Eq. (55) from exclusive B —> n£v decay as well as using the determination of \V u b\ from inclusive 
B X u lv decay [66]. In these plots, the outer error band denotes the total uncertainty in dB/dq 2 , while the inner error band 
removes the error contribution from \V u b\- 


metry breaking due to the large r mass, or more gener¬ 
ally sensitive to any Standard-Model extensions with new 
scalar currents. Moreover, the ratio of /Lt/r differential 
decay rates [67] 


and B s —> Klv form factors in Tables XI and XII. In¬ 
tegrating over the kinematically allowed ranges, we ob¬ 
tain the following Standard-Model predictions for R T Jj l = 

r ( B (s) Ptu)/Y{B^ P^v): 


v r/^ ( 2 , dT ( B (s) -t PTv)/dq 2 
p [<1 j_ dT(5 (s) -> P^)/dq 2 


(63) 


RjJ* 1 = 0.69(19), (64) 

r t ^ = 0.77(12). (65) 


provides a precise test of the Standard Model that is in¬ 
dependent of the CKM matrix element V u b \. Figure 19 
shows the predictions for the ratios of differential branch¬ 
ing fractions using our determinations of the B —>■ wiu 


The three-body final state in B ( s ) —> P£v decay also 
enables one to construct and study observables that de¬ 
pend on the kinematics of the decay products. Such 
angular observables are particularly sensitive to possi- 
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FIG. 19. Standard-Model ratio of differential branching frac¬ 
tions lZ T p(q 2 ) using our determinations of the B —> nip and 
B s —» Kip form factors in Tables XI and XII. 


ble right-handed currents. An example is the forward- 
backward difference 3 


a b (°>- 

^FB 


>Piv 


(q 2 ) 




d cos Of, 


d 2 T(B {s) P^p) 
dq 2 d cos 9( 


( 66 ) 

where 9i is the angle between the charged-lepton and 
_B( s )-meson momenta in the q 2 rest frame. In the Stan¬ 
dard Model, the forward-backward difference is given 
by [67] 


A B(ar 

^FB 


>Ptv 


(q 2 ) = 


G 2 f \V ub \ 2 
32ir 3 M Ba 



X 



M 2 K )f + (q 2 )f 0 (q 2 ) .(67) 


Integrating over the full kinematic ranges we obtain 

rQ max „ 

/ dq 2 A^^(q 2 )/\V ub \ 2 = 0.028(19) ps" 1 , (68) 

f^max 

/ dq 2 A^ T, '(q 2 )/\V ub \ 2 = 1.08(35) ps" 1 , (69) 

Jm J 

f^rnax 

/ dq 2 A%C ^{q 2 )/Kbl 2 = 0.0175(87) ps" 1 (70) 


dq 2 A^ KTV (q 2 )/\V ub \ 2 = 0.93(18) ps" 1 , (71) 


where the fi results are much smaller than the r results 
due to helicity suppression from the small muon mass. 

The normalized forward-backward asymmetry is par¬ 
ticularly interesting because it removes the ambiguity 
from \V u b\i and the theoretical form-factor uncertainties 
cancel to some degree: 




(s)- 

FB 


v 


r9 m ax J 2 A 

Jml ^ 


H*Y 

FB 


¥Ptv 


( q 2 ) 


f%r d< i 2 d P B G) p ^)/dq 2 


(72) 


The Standard-Model predictions for A^ B ^ Pl ’' (q 2 ) with 
£ = H,t are given in Fig. 21. The percentage un¬ 
certainty in A F ( B (q 2 ) is indeed smaller than in 
(q 2 )- Integrating A F ( B ^ Pei/ (q 2 ) over the full 
kinematic range we obtain 

A FB ntlv = 0.0044(13), (73) 

A F b WTV = 0.252(12), (74) 

Ayjp Kt ‘ V = 0.0039(11) (75) 

Ap^ Krv = 0.2650(79), (76) 

with errors of about 30% and 3-5% for the /i and r final 
states, respectively. 


VII. RESULTS AND CONCLUSIONS 


Figure 20 shows the Standard-Model predictions for 
(q 2 ) with £ = n, t using our determination of 
114b | from B —► nlv decay (Eq. (55)) and the inclusive 
determination of |14fc| from B —► X u £v decay [66]. Again, 
the theoretical errors are sufficiently small at large q 2 
that an experimental measurement with commensurate 
precision could distinguish between the two predictions. 


! This quantity is sometimes referred to as an “asymmetry” in the 
literature, but does not satisfy the convention that its magnitude 
is bounded by unity. 


We have calculated the form factors for B —> i t£p and 
B s —> K£u semileptonic decay in dynamical lattice QCD 
using (2+1) flavors of domain-wall light quarks and rel¬ 
ativistic b quarks. We then extended our results from 
the simulated range of lattice momenta to q 2 = 0 using 
the model-independent 2 -expansion based on analyticity 
and unitarity. We obtain our preferred results for the 
form factors f+(q 2 ) and fo(q 2 ) using the BCL form of 
the 2 -expansion [25] and imposing the kinematic con¬ 
straint /+(0) = /o(0) and a constraint on the sum of the 
coefficients for f + (q 2 ) based on heavy-quark power count¬ 
ing [62], The resulting BCL 2 -coefficients for f+(q 2 ) and 
fo(q 2 ) for B —> 7 xtv are given, along with their corre¬ 
lation matrix, in Table XI, while the BCL 2 -coefficients 
for B s —» Ktv are in Table XII. These results can be 
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FIG. 20. Standard-Model predictions for the forward-backward asymmetry for B s —> K^lu (left) and B s —> Ktu (right) using 
our determinations of the B s — > K£u form factors in Table XII. The plots show predictions for Afb using our determination of 
\V u b\ in Eq. (55) from exclusive B — > -niv decay as well as using the determination of \V u b\ from inclusive B — > X u lv decay [66]. 
In these plots, the outer error band denotes the total uncertainty in dB/dq 2 , while the inner error band removes the error 
contribution from \V u b\- 




q 2 [GeV 2 ] q 2 [GeV 2 ] 


FIG. 21. Standard-Model predictions for the normalized forward-backward asymmetry Afb for decays to muon (left) and 
r-lepton (right) final states using our form-factor determinations in Tables XI and XII. 


combined with current and future experimental measure¬ 
ments of the B — > 7 t£v and B s — > K£v branching fractions 
to obtain the CKM matrix element |I4b|- 

Figure 23, top, compares our B —» 7 tIv form-factor 
determinations with other theoretical calculations from 
light-cone sum rules (LCSR) [68, 69], NLO perturba¬ 
tive QCD (pQCD) [70], and (2+l)-flavor lattice QCD 
(LQCD) [26, 27]. Both of the earlier lattice calculations 
use staggered light quarks. The HPQCD collaboration 
uses NRQCD b quarks, while the Fermilab Lattice & 
MILC collaborations use relativistic b quarks with the 
Fermilab interpretation. Our result for the form factor 


f^ 77 agrees with earlier determinations, and is slightly 
more precise. Our scalar form factor f^ 77 is lower than 
the HPQCD result (although by less than 2cr), but we 
note that their calculation used only a single lattice spac¬ 
ing. We also agree with the most recent light-cone-sum- 
rule prediction for the B —> 7 xtv form factor from Imsong 
et al. [69], who present the first extrapolation of LCSR 
results from low q 2 to g^ ax using the z expansion. 

We fit our results for the B —► niv form factors to¬ 
gether with the experimentally measured decay rates 
from BaBar [2, 4] and Belle [3, 5], leaving the relative nor¬ 
malization as a free parameter, to determine the CKM 
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HFAG inclusive 
FLAG (N f = 2+1) 
CKMfitter Group 
UTfit Collaboration 

HPQCD 2006 (q 2 > 16GeV 2 ) 
FNAL/MILC 2009 (BCL z-fit) 
This work 


FIG. 22. Determinations of \V u b\ from Table XIV. For points 
with double error bars, the inner error bars are experimental 
while the outer error bars show the total experimental plus 
theoretical uncertainty added in quadrature. 
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CKM unitarity 
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- \ —B—I B —' niv — 

I->-1 

■ i I i i i i I i i i i I i i i i I i i i i I 

3.0 3.5 4.0 4.5 5.0 

IvJ x 10 3 


matrix element \V u b\. We obtain 

\V ub \ = 3.61(32) x lO" 3 , (77) 

where the error includes both theoretical and experimen¬ 
tal uncertainties. Table XIV and Figure 22 compare the 
determination of \ V ub \ using our lattice form-factor result 
with determinations using other theoretical calculations 
of the B -A niv form factor, as well as with determina¬ 
tions from inclusive B —> X u iv decay, B -A tv leptonic 
decay, and predictions from CKM unitarity. Our \V ub \ 
result agrees with other lattice determinations, as well 
as with the less precise determination from B -A tv de¬ 
cay, and with the more precise predictions from CKM 
unitarity. The central value is higher, however, than the 
one obtained from the FNAL/MILC calculation so the 
tension between our result and the determination from 
inclusive B -A X u iv decay is less, only about 2er. 

Figure 23, bottom, compares our B s —> Kiv form- 
factor determinations with the theoretical calculation 
by the HPQCD collaboration using (2+l)-flavor lat¬ 
tice QCD with staggered light quarks and NRQCD b 
quarks [30], as well as with predictions at q 2 = 0 from 
the light-cone sum rules [73], NLO perturbative QCD [74] 
and the relativistic quark model (RQM) [75]. The lattice 
results agree in the range of simulated lattice data, and 
have similar precision, but diverge slightly when extrapo¬ 
lated to lower q 2 . Even at q 2 = 0, however, the predicted 
form factors f+“ K { 0) differ by only 1.9cr. 

The semileptonic decay B s -A Kiv has not yet been 
observed experimentally. We can therefore make predic¬ 
tions for Standard-Model observables for the decay pro¬ 
cesses B s —k K^v and B s —k Ktv. Using our results 
for the B s -A Kiv form factors in Table XII and our de¬ 
termination of \V ub \ from B —► nil/ given above, we cal¬ 
culate the differential branching fractions and forward- 
backward asymmetries. Our lattice form-factor deter¬ 
minations at q 2 > 13GeV 2 are sufficiently precise that 
future experimental measurements of the B s -A Kiv dif¬ 
ferential branching fraction in this range with similar un¬ 


certainties will be able to distinguish between Standard- 
Model predictions using \V ub \ from inclusive B —k X u iv 
and exclusive B —> niv semileptonic decays, and thereby 
weigh in on the current ~ 3cr disagreement between the 
two determinations. 

We also calculate the ratio of /x-to-r differential de¬ 
cay rates, and the normalized forward-backward asym¬ 
metries, for both B —> niv and B s -A Kiv. Because 
these quantities are independent of \V ub \, they potentially 
provide more stringent tests of Standard-Model exten¬ 
sions such as ones that give rise to new scalar or right- 
handed currents. In practice, it will likely be difficult for 
LHCb to measure B s decays with r leptons in the final 
state, so the B s -a K^v predictions in Sec. VIB may be 
most useful for the foreseeable future. Belle II, however, 
should observe the decay B -A niv with both i = //, r 
final states, and we anticipate that they will measure the 
forward-backward asymmetries and eventually the /x/r 
ratio for this decay. Future measurements of these ob¬ 
servables will be especially important given the current 
> 3cr discrepancy observed in R{D) and R(D*) for the 
similar decays B -A D^iv [76, 77]. 

Our results for the B -A niv and B s —> Kiv pro¬ 
vide important independent checks of existing calcula¬ 
tions using staggered light quarks. Such confirmation 
is especially important given the present approximately 
3cr tension between V ub obtained from inclusive and 
exclusive semileptonic B decays. Currently the preci¬ 
sion of our determination of the form factors is limited 
by statistics and by the relatively large chiral extrapo¬ 
lation to the physical light-quark mass. To address the 
chiral-extrapolation error, we are presently analyzing the 
RBC/UKQCD Mobius domain-wall + Iwasaki ensem¬ 
ble [78-80] with a lattice spacing close to the coarser 
value a ~ 0.11 fm used in our current analysis, but 
with M n « 140 MeV. We are also using all-mode av¬ 
eraging [81, 82] to reduce the statistical errors on the 
individual numerical data points, and expect some reduc¬ 
tion in the statistical errors. With these improvements 
we anticipate a reduction in the current form-factor er¬ 
rors. Further, our future physical-mass results will also 
include a determination of the tensor form factor, which 
will enable a calculation of the Standard-Model rate for 
the rare decay B -A ni + i~ and similar processes. 

Because we present our results for the B -A niv and 
B s -A Kiv form factor as coefficients of the BCL z- 
parametrization and the matrix of correlations between 
them, our form factors can be combined with new ex¬ 
perimental measurements (and even with other lattice 
form-factor calculations) to further improve V ub in the 
future. In particular, our B s —> Kiv form-factor results 
will enable an alternate determination of \V ub \ once the 
process has been observed in experiment. More generally, 
our form-factor results in Tables XI and XII can be used 
to compute all possible Standard-Model observables for 
these decays whenever they are needed for comparison 
with experiment. 
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TABLE XIV. Determinations of \V u b\- Top panel: results from inclusive B —» X u lv decay [66] and B —» tv leptonic decay [11]. 
Middle panel: predictions from CKM unitarity [71, 72]. Bottom panel: results from exclusive B —» nip decay using form factors 
from (2+l)-flavor lattice QCD [26, 27, 66]. Errors shown are either the total uncertainty or the experimental and theoretical 
uncertainties, respectively. 



From 

\Vub\ x 1(U 

HFAG inclusive average [56] 
FLAG (Nf = 2 + 1) [11] 

B -+ X u lu 

B —¥ TV 

4.40(15)(20) 

4.18(52)(9) 

CKMfitter Group [71] 

UTfit collaboration [72] 

CKM unitarity 3.435(tg4°) 
CKM unitarity 3.63(12) 

HPQCD (HFAG q 2 > 16GeV 
FNAL/MILC (HFAG BCL z- 
This work 

2 ) [26, 66] B -+ n£p 
-fit) [27, 66] B -+ nhs 

B —> 7 t£v 

3-52(8) (Zt®o) 

3.28(29) 

3.61(32) 
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FIG. 23. Top: theoretical calculations of the B —> niv form factors from light-cone sum rules [68, 69], NLO perturbative 
QCD [70], and (2+l)-flavor lattice QCD [26, 27]. Bottom: theoretical calculations of the B a —¥ K£v form factors from QCD 
models [73-75] and (2+1)-flavor lattice QCD [30]. In all plots, the predictions for f+(q 2 = 0) are displayed with a slight 
horizontal offset for clarity. 
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Appendix A: Correlator fit results 

Here we summarize the results for pion, kaon, B and 
B s meson masses and three-point ratios from the corre¬ 
lator fits described in Secs. Ill A and IIIB. 


Appendix B: 2-fit results 

Here we present the complete results for the 2 -fits de¬ 
scribed in Secs. V B and VI A. Tables XIX and XX show 
the results of fits to only our lattice form factors, while 
Tables XXI and XXII show the result of joint fits to the 
lattice and experimental data. 
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TABLE XV. Pion and kaon masses on all ensembles. See Sec. Ill A for details. 


mi 

[^minj ^max] 

Pion 

aMn 

X 2 /dof 

P 

[^min? ^max] 

Kaon 

clMk 

X 2 /dof 

P 

0.005 

[12:23] 

0.18959(53) 

0.83 

61% 

[12:23] 

0.31287(45) 

1.31 

21% 

0.01 

[12:23] 

0.24305(48) 

0.76 

68% 

[12:23] 

0.33352(44) 

0.82 

62% 

0.004 

[16:30] 

0.12611(51) 

1.08 

37% 

[16:30] 

0.23249(42) 

1.35 

17% 

0.006 

[16:30] 

0.15207(36) 

0.41 

97% 

[16:30] 

0.24189(33) 

0.58 

88% 

0.008 

[16:30] 

0.17265(42) 

1.05 

40% 

[16:30] 

0.24854(39) 

1.06 

39% 


TABLE XVI. B- and B s -meson masses on all ensembles. See Sec. Ill A for details. 


mi 

[^minj ^max] 

B meson 
aMg 

X 2 /dof 

P 

[^•min? ^max] 

B s meson 
aMg 

X 2 /dof 

P 

0.005 

[7:30] 

3.0638(13) 

0.92 

57% 

[10:30] 

3.1020(12) 

0.46 

98% 

0.01 

[7:30] 

3.0727(13) 

0.74 

80% 

[10:30] 

3.1028(13) 

1.13 

31% 

0.004 

[9:30] 

2.3203(14) 

1.40 

11% 

[13:30] 

2.3509(11) 

1.53 

7% 

0.006 

[9:30] 

2.3241(10) 

0.78 

74% 

[13:30] 

2.3520(09) 

0.57 

92% 

0.008 

[9:30] 

2.3274(12) 

0.85 

66% 

[13:30] 

2.3533(12) 

0.95 

51% 


TABLE XVII. Three-point correlator ratios R 3 J,. on all ensembles. Results are shown for pion momenta = ( 2-ivn/L ) 2 
through n 2 = 3. 


mi 

n 2 

[^min? ^max] 

Mg^Rfj/pi 

X 2 /dof 

P 

[^min? ^max] 

Mg] /2 R$£ 

X 2 /dof 

P 

0.005 

0 





[6:10] 

0.2523(36) 

1.49 

20% 


1 

[6:10] 

1.086(31) 

1.19 

31% 

[6:10] 

0.2034(52) 

0.16 

96% 


2 

[6:10] 

0.827(46) 

1.71 

15% 

[6:10] 

0.1819(93) 

0.38 

82% 


3 

[6:10] 

0.755(78) 

0.60 

66% 

[6:10] 

0.154(17) 

0.62 

65% 

0.01 

0 





[6:10] 

0.2513(35) 

2.03 

9% 


1 

[6:10] 

1.049(23) 

1.07 

37% 

[6:10] 

0.2060(43) 

1.17 

32% 


2 

[6:10] 

0.798(30) 

0.53 

72% 

[6:10] 

0.1837(73) 

0.55 

70% 


3 

[6:10] 

0.706(47) 

0.88 

47% 

[6:10] 

0.180(13) 

0.04 

100% 

0.004 

0 





[8:13] 

0.3679(79) 

1.25 

28% 


1 

[8:13] 

1.546(71) 

0.88 

50% 

[8:13] 

0.288(12) 

1.58 

16% 


2 

[8:13] 

1.100(92) 

0.40 

85% 

[8:13] 

0.256(18) 

1.93 

9% 


3 

[8:13] 

0.96(20) 

0.33 

89% 

[8:13] 

0.186(32) 

0.70 

63% 

0.006 

0 





[8:13] 

0.3528(51) 

1.19 

31% 


1 

[8:13] 

1.529(42) 

1.01 

41% 

[8:13] 

0.2739(63) 

0.30 

91% 


2 

[8:13] 

1.219(58) 

0.41 

84% 

[8:13] 

0.240(11) 

0.76 

58% 


3 

[8:13] 

0.921(91) 

0.46 

80% 

[8:13] 

0.224(21) 

1.14 

34% 

0.008 

0 





[8:13] 

0.3432(58) 

1.45 

20% 


1 

[8:13] 

1.483(38) 

1.43 

21% 

[8:13] 

0.2829(67) 

1.02 

40% 


2 

[8:13] 

1.183(53) 

0.63 

68% 

[8:13] 

0.2349(98) 

1.01 

41% 


3 

[8:13] 

0.888(89) 

0.52 

76% 

[8:13] 

0.211(18) 

0.20 

96% 
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TABLE XVIII. Three-point correlator ratios on all ensembles. Results are shown for kaon momenta p% = (27rn/L) 2 

through n 2 = 4. 


mi 

n 2 

[^min? ^max] 

K^Kf/PK 

X 2 /dof 

P 

[^minj ^max] 

M^Rlf 

X 2 /dof 

P 

0.005 

0 





[6:10] 

0.2394(27) 

0.99 

41% 


1 

[6:10] 

0.984(16) 

2.17 

7% 

[6:10] 

0.2036(27) 

0.60 

66% 


2 

[6:10] 

0.763(16) 

2.53 

4% 

[6:10] 

0.1791(36) 

0.46 

77% 


3 

[6:10] 

0.623(22) 

1.60 

17% 

[6:10] 

0.1609(57) 

0.52 

72% 


4 

[6:10] 

0.543(31) 

0.98 

42% 

[6:10] 

0.1568(88) 

0.68 

61% 

0.01 

0 





[6:10] 

0.2421(28) 

2.35 

5% 


1 

[6:10] 

0.963(15) 

0.95 

44% 

[6:10] 

0.2060(30) 

1.00 

40% 


2 

[6:10] 

0.749(16) 

0.13 

97% 

[6:10] 

0.1847(40) 

0.17 

95% 


3 

[6:10] 

0.634(21) 

0.38 

83% 

[6:10] 

0.1737(61) 

0.16 

96% 


4 

[6:10] 

0.573(30) 

0.30 

88% 

[6:10] 

0.1664(81) 

0.12 

98% 

0.004 

0 





[8:13] 

0.3264(49) 

1.14 

34% 


1 

[8:13] 

1.340(26) 

0.57 

72% 

[8:13] 

0.2751(51) 

1.25 

28% 


2 

[8:13] 

1.039(26) 

1.14 

33% 

[8:13] 

0.2444(64) 

1.84 

10% 


3 

[8:13] 

0.830(39) 

1.65 

14% 

[8:13] 

0.2148(89) 

1.80 

11% 


4 

[8:13] 

0.690(50) 

0.88 

49% 

[8:13] 

0.195(13) 

1.18 

32% 

0.006 

0 





[8:13] 

0.3312(35) 

1.12 

34% 


1 

[8:13] 

1.336(20) 

1.02 

40% 

[8:13] 

0.2770(36) 

1.19 

31% 


2 

[8:13] 

1.068(20) 

0.71 

61% 

[8:13] 

0.2419(46) 

0.60 

70% 


3 

[8:13] 

0.868(26) 

1.49 

19% 

[8:13] 

0.2168(68) 

0.53 

75% 


4 

[8:13] 

0.737(37) 

0.57 

72% 

[8:13] 

0.206(10) 

0.37 

87% 

0.008 

0 





[8:13] 

0.3323(47) 

1.28 

27% 


1 

[8:13] 

1.332(25) 

1.21 

30% 

[8:13] 

0.2848(47) 

1.16 

32% 


2 

[8:13] 

1.056(25) 

0.86 

51% 

[8:13] 

0.2476(58) 

1.23 

29% 


3 

[8:13] 

0.856(34) 

1.17 

32% 

[8:13] 

0.2264(87) 

0.80 

55% 


4 

[8:13] 

0.736(45) 

0.83 

52% 

[8:13] 

0.200(13) 

0.61 

69% 
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TABLE XXI. Results for K = 3 2 -fits of the B —> nIv experimental data to the BCL parametrization, Eqs. (44) and (45). 



IKil&f x 10 3 



\V u b\ 2 J2 B rnnb m b n X 10® 

X 2 /dof 

V 

All 

1.517(32) 

-1.03(19) 

-1.18(61) 

0.23(12) 

1.36 

5% 

BaBar 2010 [2] 

1.361(75) 

-0.91(47) 

0.5(1.5) 

0.06(11) 

1.98 

11% 

BaBar 2012 [4] 

1.497(58) 

-0.33(44) 

-3.5(1.3) 

1.19(77) 

0.45 

91% 

Belle 2010 [3] 

1.602(61) 

-1.31(27) 

-0.8(0.9) 

0.21(12) 

1.19 

29% 

Belle 2013 [5] 

1.562(84) 

-1.88(50) 

1.3(1.6) 

0.40(53) 

1.24 

22% 


TABLE XXII. Determinations of \V u b\ from combined fits of B —> ir£v experimental and lattice form-factor data to the BCL 
parametrization, Eq. (44), for different truncations K. 



K 

|14b| x 10 3 

6<°> 

b?/b^ 

b^/bf 

6l 3) /&l 0) 

bf/b { ? 

) y Bmnbmbn 

X 2 /dof 

V 

All 

3 

3.61(32) 

0.422(35) 

-1.07(17) 

-1.08(56) 



0.0162(86) 

1.32 

7% 


4 

3.61(32) 

0.421(37) 

-1.07(18) 

-1.0(1.6) 

-0.6(4.1) 


0.015(19) 

1.35 

5% 


5 

3.61(32) 

0.421(37) 

-1.10(27) 

-1.0(1.6) 

-1(12) 

-5(26) 

0.2(2.4) 

1.37 

4% 

BaBar 2010 [2] 

3 

3.31(36) 

0.414(37) 

-1.03(31) 

0.9(1.1) 



0.009(11) 

1.35 

24% 


4 

3.29(35) 

0.428(41) 

-1.04(30) 

-1.4(2.6) 

-7.5(7.2) 


0.47(97) 

1.48 

20% 


5 

3.29(35) 

0.428(41) 

-1.21(39) 

-1.0(2.7) 

17(15) 

-20(39) 

7(20) 

1.82 

14% 

BaBar 2012 [4] 

3 

3.52(35) 

0.436(38) 

-0.79(28) 

-2.20(86) 



0.043(31) 

0.61 

82% 


4 

3.56(36) 

0.425(40) 

-0.81(29) 

-0.5(2.4) 

-6.2(6.5) 


0.24(54) 

0.61 

81% 


5 

3.54(36) 

0.428(40) 

-1.04(36) 

-0.3(2.4) 

9(15) 

-41(35) 

15(27) 

0.54 

85% 

Belle 2010 [3] 

3 

3.90(37) 

0.411(35) 

-1.29(23) 

-0.86(79) 



0.0141(82) 

1.01 

43% 


4 

3.93(38) 

0.403(36) 

-1.27(23) 

0.6(2.0) 

-4.8(6.0) 


0.14(38) 

1.05 

40% 


5 

3.94(38) 

0.401(36) 

-1.41(34) 

1.0(2.2) 

2(13) 

-24(36) 

4(13) 

1.12 

34% 

Belle 2013 [5] 

3 

3.88(44) 

0.397(37) 

-1.52(36) 

0.2(1.2) 



0.0109(56) 

1.17 

27% 


4 

3.84(42) 

0.425(41) 

-1.59(32) 

-5.0(2.6) 

-17.0(7.6) 


2.6(2.5) 

1.01 

45% 


5 

3.86(42) 

0.429(41) 

-1.28(38) 

-6.2(2.7) 

-2(15) 

73(41) 

41(52) 

0.95 

51% 
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